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Abstract 

We consider a nonlinear Schrodinger equation with a bounded local poten- 
■ tial in K 3 . The linear Hamiltonian is assumed to have three or more bound 

states with the eigenvalues satisfying some resonance conditions. Suppose that 
the initial data is localized and small of order n in H , and that its ground state 
\q • component is larger than w?~ e with e > small. We prove that the solution 

will converge locally to a nonlinear ground state as the time tends to infinity. 

o 

(N : 1 Introduction 

o 



X 



Consider the nonlinear Schrodinger equation 



-g: z<9 t V> = (-A + V)V + Ah/>| 2 V, ^(t = 0)=Vo, (1.1) 

a 

where V = V(x) is a smooth localized real potential, A = ±1 and ip = ip(t,x) : 
Rxl 3 — > C is a wave function. For any solution ip(t) e H 1 (M. 3 ) the L 2 -norm and 
the Hamiltonian 



WM = f\w? + \v\^ + \\Wdx (1.2) 

are constant for all t. The global well-posedness for small solutions in if 1 (M 3 ) can be 
proven using these conserved quantities and a continuity argument, no matter what 
the sign of A is. We assume that the linear Hamiltonian H := — A+V has iV+1 simple 
eigenvalues e < e\ < ■ ■ ■ < with normalized eigenvectors (fik, k — 0, 1, . . . , N, 
where N > 2. These eigenvalues are assumed to satisfy some resonance conditions to 



be specified later on. The nonlinear bound states to the Schrodinger equation (|P 
are solutions to the equation 

(-A + V)Q + \\Q\ 2 Q = EQ, (1.3) 
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for some E. They are critical points to the Hamiltonian Tl[i/)] defined in (|1.2j) subject 
to the constraint of fixed L 2 -norm. For any nonlinear bound state Q = Qe{x), 
the function xjj(t, x) = Q(x)e~ tEt is an exact solution to the nonlinear Schrodinger 
equation fll.ip . We may obtain A^ + 1 families of such nonlinear bound states by 
standard bifurcation theory, corresponding to the N + 1 eigenvalues of the linear 
Hamiltonian. For any E sufficiently close to eo so that E — eo and A have the same 
sign, there is a unique positive solution Q = Q E to ( |1.3|) which decays exponentially 
as x — > oo. We call this family the nonlinear ground states and we refer to it as 
{Qe}e- Similarly, for each k > there is a nonlinear excited state family {Qk,E k } Ek 
for Ek near e^. These solutions are small and ||Qfc,£ fe || ~ \Ek — e^ 1 / 2 . See Lemma \2,.1[ 
Our goal is to understand the long-time dynamics of the solutions at the presence 
of nonlinear bound states. One first considers the stability of nonlinear ground states. 
There are two different concepts: orbital stability and asymptotic stability. It is well- 
known that nonlinear ground states are orbitally stable in the sense that the difference 



inf \\m-Qse 



iO I 



\L 2 



remains small for all time t if it is initially small. On the other hand, one expects that 
the difference actually approaches zero locally, as the majority of the difference is a 
dispersive wave which escapes to infinity. Hence one expects that it is asymptotically 
stable in the sense that 

\\m-QE(t)e te{t) \\ L2 ^0 

loc 

as t — > oo, for a suitable choice of E(t) and O(t). Here ||-|| L 2 denotes a local L 2 

loc 

norm, a precise choice will be made later on in ( |1.8[ ). One also wants to determine 
the decay rate and whether E(t) has a limit. A more difficult problem, which one 
studies next, is the asymptotic dynamics of the solutions when the initial data are 
away from nonlinear ground states. As in the previous problem, one wants to identify 
their local behavior. 

If —A + V has only one bound state, i.e., with no excited states, the asymptotic 
stability is proved in |16|]. The solution eventually settles down to some ground state 
Qeoo with E^ close to E and the local difference is bounded by Ct -3 / 2 , as the decay 
rate of the free evolution e ltA . This result is extended in [10] to all small initial data, 
not necessarily near ground states. See also |l|, [|. 



Suppose —A + V has two bound states. We proved in |T9| that the evolution 
with initial data tpo near some Qe will eventually settle down to some ground state 
Qeov with Eoo close to E. The local difference is, however, only bounded by Ct^ 1 ^ 2 
due to the persistence of the excited state. The key mechanism here is the resonance 
decay of the excited state due to resonance with the continuous spectrum. See P, [| 
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for an one dimensional equation, || its extension to higher dimensions, and W7\ for 
real-valued nonlinear Klein-Gordon equations. 

The problem becomes more delicate when the initial data are away from ground 
states. Based on physical intuition, one expects that all solutions near excited states 
decay to ground states unless initially they are exactly excited states. However, we 
proved that there exists a family of "finite co-dimensional manifolds" in the space 
of initial data so that the dynamics asymptotically converge to some excited states. 
Outside a small neighborhood of these manifolds, the asymptotic profiles are given 



by some ground states f20 |. We further showed [22] that there are exactly three 
asymptotic profiles: vacuum, excited states or ground states. The last problem is 
also considered in |18| . Earlier works concerning related linear analysis were obtained 
in l||i§. 



In this paper, we extend the results in [|T^, |20[ to the case when — A + V has three 
or more bound states. 

Our assumptions on the operator —A + V are as follows: 
Assumption AO: —A + V acting on L 2 (M 3 ) has N + 1 simple eigenvalues eo < e\ < 
. . . < e^v < 0, N > 2, with normalized eigenvectors (f)Q, . . . ,4>n- 

Assumption Al: V(x) is a real-valued function. For XQ 2 E sufficiently small, the 
bottom of the continuous spectrum to —A + V + XQ%, 0, is not a generalized eigen- 
value, i.e., not an eigenvalue nor a resonance. Also, we assume that V satisfies 
the assumption in p4j so that the W k,p estimates k < 2 for the wave operator 
Wh = lim^oo e LtH °e ttA hold for k < 2, i.e., there is a small s± > such that, 

\V a V(x)\ < C(x)' 5 ' Sl , for |a| < 2. 

Also, the functions (x • V) k V, for k = 0, 1,2,3, are —A bounded with a —A bound 
less than 1: 

|| (x • V) k V<p\\ 2 < s 2 ||-A0|| 2 + C U\\ 2 , s 2 < 1, k = 0, 1, 2, 3. 
Assumption A2: Resonance condition. Denote e^o = &k ~ We assume that 

eo<2ei. (1.4) 
Hence 2eko > l e o| f° r all > 0. We further assume that, for some small Sq > 0, 

7o= inf lim Lim ( (f> <f>l, — — — — — — ;P c o 0l)>O. (1.5) 

i<k<N r^o+ \ —A + V + e — 2et — s — ri J 

s|<s 

Assumption A3: No-resonance condition. Let j max = 3. For all j = 2, . . . ,j max and 
for all kx, . . . ,kj,li, . . . , lj G {0, . . . , iV} with {k\, . . . , kj} ^ {1%, . . . , lj} as sets with 
multiplicities, (e.g., {0, 0, 1} ^ {0, 1, 1}), 

ejfcj H h e k , ^e h -\ h e/.. (1.6) 
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Assumption Al contains some standard conditions to assure that most tools for 
linear Schrodinger operators apply. These conditions are certainly not optimal. (For 
example, it is sufficient to assume is not a resonance or eigenvalue of Hq, which 
implies the same statement for -A + V + XQ 2 . See plj.) The main assumption in A2 



is the condition e < 2e 1 . Since the expression for j is quadratic, it is non- negative 
and 70 > holds generically. The condition eo < 2e± states that the energies of the 
excited states are closer to the continuum spectrum than to the ground state energy. 
It guarantees that, for each k > 0, 1e^ — eo > becomes a resonance in the continuum 
spectrum of Hq. (Hq + eo — 2et is not invertible in L 2 .) This resonance produces the 
main relaxation/growth mechanism. If this condition fails, the resonance occurs in 
higher order terms and a proof of relaxation will be much more complicated. Also, 
the rate of decay will be different. Assumption A3 is a new condition to avoid direct 
resonance between the eigenvalues. It is trivial if AT = 0, 1. It holds true generically 
and is often seen in dynamical systems of ODE's. See Example 2 at the end of this 
section for what may happen if this assumption fails. If we relax the assumption 
eo < 2ei, we may need to increase j max - 

Denote by L 2 the weighted L 2 spaces (r may be positive or negative), 

L 2 (R 3 ) = {0 G L 2 (R 3 ) : (1 + x 2 ) r/2 <J) G L 2 (R 3 )} . (1.7) 



Fix ri > 3 large enough, to be determined by estimates (|2]5|) and (|2.30|) . We denote 
by Lf oc , p = 1, 2, the local LP spaces given by the norm 



loc 



l/p 

(1 + \x\)- 2ri \4>\ p d 3 x\ . (l.i 



Now we are ready to state our main theorem. 

Theorem 1.1 Suppose the assumptions A0-A3 on Hq = —A + V hold and let e > 
be any small constant. There is a small constant no > such that the following holds. 
Let ijj(t,x) be a solution of ( |1 . 1| ) with initial data ipo decomposed as 

= Xq4> + x° x (j)i H h x° N 4> N + £ 

with respect to H , where x° k G C, k — 0, 1, . . . ,N, and £ G H C (LL ). Suppose that 

||^o||y =n, 0<n<n , \x° \ > n 3 " e , ||£ ||y < n/2, (1.9) 

where Y is the space for initial data, 

Y = H l ni^R 3 ). (1.10) 
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Then there exists an with \[Qe x \[i? ~ n an d a real function 6(t) = —E^t + 
O(logt) such that 

Wm-QE^w^ <c 2 {i+t)- i i\ (1.11) 

loc 

/or some constant C2 > depending on n. Suppose, furthermore, maxf =1 > 
n/100, we a/so /iave a /ower bound 

Wm-QE^w^ ^c^i+t)- 1 ' 2 , (1.12) 

loc 

/or some constant C\ > depending on n. 

The condition > n 3_e is certainly not optimal and can be greatly relaxed. 
It ensures that -?/>o is away from nonlinear excited states and that the dispersion 
component, £05 is n °t extremely large compared with x®. This condition, however, 
still allows the ground state component to be much smaller than other components 
and exhibit the main phenomena. 

Depending on the relative sizes of the bound states, there are three regimes: 

I. when an excited state is dominant; 

II. when the ground state and the excited states are comparable; 

III. when the ground state is dominant. 

Because the dominant terms are different in different regimes, the natural linear 
operators and the corresponding decompositions of the wave function are different. 
In regime II we can use Ho as the linear operator and decompose a function ip e L 2 

as 

^ = aj o o H Vx N (t) N + i, (1.13) 

with Xk £ C and £ G H c (Hq). When the function if) is close to a nonlinear ground 
state Qe, i.e., in regime III, it is natural to use C = Ce, the linearized operator 
around Qe, 

Ch = -i{(-A + V - E + 2XQ 2 )h + XQ 2 h} , (1.14) 

and decompose the wave function according to the spectral decomposition with re- 
spect to C Notice that C is not self-adjoint. Similarly, we can use linearized operators 
around excited states in regime I. 

The picture of the dynamics is as follows. Suppose the ground state component is 
initially of order n 3 ~ e and the solution is in regime II. The ground state component 
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will gradually grow while the excited states gradually decay, until the solution enters 
regime III, i.e., when the ground state component becomes much larger than other 
components. This time interval is called the transition regime. After entering regime 
III, the solution will converge locally to some nonlinear ground state as time tends to 
infinity, with the excited state components vanishing and the dispersion component 
escaping to infinity. This time interval is called the stabilization regime. 

As indicated above, we will use different operators and coordinate systems for these 
two regimes. Besides technical problems associated with changing coordinate systems, 
there is an intrinsic difficulty related to the time reversibility of the Schrodinger 
equations; It is not sufficient to control only the usual Sobolev space norms at the 
time of coordinates changing. To overcome this difficulty, we will use a concept 
introduced in |2(| p2| , the out-going estimates, to capture the time-direction sensitive 
information of the dispersive waves. 

As it will be seen from examples below, there are energy transfers from higher 
modes to lower modes. Although all excited states eventually vanish, in an interme- 
diate time an excited state may actually grow because it gains more energy than it 
loses. This complicates the analysis: one cannot prove the decay of each individual 
excited state for all time. Instead, we have derived some monotonicity formulas for 
their sum, see (|4.57|) . Because the energy transfer between excited states is relatively 



small in the stabilization regime, the phenomenon mentioned above is only apparent 
in the transition regime. 

We now give two examples illustrating some phenomena of many bound states, 
in particular the relaxation/growth mechanism. We first recall the concept of phase 
factor. In our analysis it is essential to identify the main oscillation factors of various 
terms. If we decompose the solution ip(t) according to ( |1.13| ), the equation for a 
component Xk(t) is 

ix k = e k x k + (4> k , \ip 2 4>). (1.15) 

From the linear part we find that Xk{t) has an oscillation factor e~ t£kt . We write 
Xk{t) = e~ lSkt Uk{t) and say that its phase factor is — e^. We will talk about phase 
factors of polynomials in Xk in a similar way. For example, the phase factor of x\x m Xj 
is e\ e m -\- Gj. 

Example 1. Three bound states case in transition regime 

Suppose H has three bound states and denote x{t) = xo(t), y(t) = xi(t) and 
z(t) = X2{t). The leading terms of £ are generated by bound states and are cubic in 
Xj. Those with resonant coefficient functions, i.e., those with negative phase factors, 
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give the relevant part of £, 



£, = (y 2 + z 2 + yz)x + z 2 y H (1.16) 

Here we have ignored their coefficient functions, which are non-local and complex 
because of the resonance with the continuous spectrum. Denote x — 



N 



Xj<pj 



x<fio + ytfii + z<p2- The main terms of ( |1.15p for x k are 



ix k = e k x k + (0 fc , \x 2 £ + 2\xx0 H , (1-17) 

where • • • denotes irrelevant terms. Substituting ( |1.16|) into ( |1.17|) , we get 



x= [7 2 ° 2 k| 4 + 47? 2 M 2 k| 2 + 7? 1 M 4 ] x + ... 

y= H2\z\*-H2\x\ 2 \z\ 2 -Hi\x\ 2 \y\ 2 } y + --- (lis) 

z= [-2 72 2 |x| 2 |z| 2 - 47 1 2 |x| 2 |y| 2 - 2 72 1 2 | 2 /| 2 | 2 | 2 ] z + 

Here ■ ■ • denotes irrelevant and error terms; 7^, 7° 2 , 7 22 an d 7 22 are non- negative 
constants to be defined later in (|4.26|) . These constants (except 7 22 , depending on 
whether e\ — 2e 2 < 0) are generically positive. They are the nonlinear analogue of 
the Fermi golden rule, extensively studied in, e.g., |T^, [T7|, 0, |3|, [5|, ^TJ. The irrelevant 



terms have two kinds. The first kind consists of terms with different phase factors. 
They have few effect averaging over time and can be removed using a normal form 
procedure, see Lemma |4.2|. The second kind consists of terms with same phase factors 



but with purely imaginary coefficients. They only contribute to the phase of x k , not 
to the magnitude. To illustrate further, let us assume 7^ = 7° 2 = 7 22 = 7 22 = 1 and 
denote A = \x\ 2 , B = \y\ 2 , C = \z\ 2 . By ( |1.18| ) we have 

A/2 = AABC + AB 2 + AC 2 + ■■■ 

B/2 = -AABC -2AB 2 + BC 2 + ■■■ (1.19) 
C/2 = -AABC - 2 AC 2 - 2BC 2 + ■■■ 

where • • • denotes irrelevant terms. Although the above system is accurate only in 
the transition regime, corresponding equations for other regimes are similar. From 
( [LIS ) we can read the energy transfers from higher modes to lower modes. Only half 



of the energy decrease in a higher mode goes to lower modes, while the other half goes 
to dispersion (radiation). The lowest mode (ground state) is only receiving energy 
while the highest mode (x^) is only losing energy. The intermediate modes receive 
energy from higher modes and release energy to lower modes and radiation. 

Richer phenomena occur when the constants 7 ; J m have different sizes. For example, 
suppose 7 22 is much larger than other constants and initially the energy is concen- 
trated in the third mode z. Then the second mode y will first grow exponentially, 
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acquiring energy from z, and then gradually decay Another phenomenon to be no- 
ticed is the interaction between three modes, the ABC terms in ( |1.19| ). They have 
the same coefficients ±47° 2 in ( 1.18| ). Hence this interaction is more apparent when 



7° 2 is much larger than other coefficients. This phenomenon is present only for many 
bound states case, N > 2. 



Example 2. No-resonance assumption A3 violated 

Still assume eo < 2ei and three bound states. The only possibility for assumption 
A3 to fail is 

e + 2e 2 = 3ei. (1.20) 
An example is e = —10, e 1 = —4 and e 2 = —1. Note that 

ex -2e 2 = e - 2ei < 0. (1.21) 



Hence i s generically positive. Because of (|1.20Q , when we substitute (|1.16| ) into 
( 1.17 ) we get new resonant terms in ( |1.18j ): 



x = [• • • ]x + y (z 2 y) + 



y=[...}y + z 2 (y 2 x) + xy(yz 2 ) + ■■■ (1.22) 
z= [■ • - \z + yz{xy 2 ) H 

The first group of terms on the right side are those in ( |1.18| ). The cubic terms 
in the parentheses () are from (|1.16|). For example, y 2 (z 2 y) has the phase factor 



— 2ei — (— 2e 2 + ei) = — e , the same as x, due to ( |1.20j ). Although these new terms 
on the right side have the same phase factors as the left side, their phases are not 
exactly the same. Moreover, their coefficients are not quadratic and it is unclear how 
to determine their signs. Hence it is difficult to predict the dynamics of this system. 



2 Preliminaries 

We first fix the notation. Let H k denote the Sobolev spaces W^'^IR 3 ). The weighted 
Sobolev space L 2 (M. 3 ) is defined in (p..7|). Denote by conj the conjugation operator. 
The L 2 inner product ( , ) is 

(f,g)= ! fg d 3 x. (2.1) 

For a function G L 2 , denote by (j) 1 - the L 2 -subspace {g e L 2 : (0, g) = 0}. 

In what follows we collect some facts about nonlinear bound states and linear 
analysis. Their proofs can be found in [|Tj], [2D], |2T|. Although the proofs there are for 
two bound states case, the proofs for the many bound states case are the same. 
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2.1 Nonlinear bound states 

Recall that nonlinear bound states of the equation (|1.1| ) are solutions of ( |1.3|) . They 
are critical points of the energy functional l~C[ip] defined in (|1.2|) , subject to the con- 
straint of fixed L 2 -norm. For each such solution Qe, the function ip(t, x) = Q E{x)e~ lEt 
is an exact solution of ( |1 . 1|) . Since we have N + 1 simple eigenvalues, we have N + 1 
families of corresponding nonlinear bound states. The existence and basic properties 
of these nonlinear bound states are summarized in the following lemma. They are 



proven in 19, 20| using a contraction mapping argument. 



Lemma 2.1 Suppose —A + V satisfies assumptions A0-A1. Let n$ be sufficiently 
small. For each eigenvalue e k with normalized bound state <f> k , k — 0, 1, . . . , N, there 
is a family of nonlinear bound states {Qk,E k } Ek t° (OD f or between e k and e k + Xn^ 
such that Qk,E k are real, localized, smooth, and X~ 1 (E k — e k ) > 0, 

Qk,E k = n<p k + h, h±(/) k , h = 0(n 3 )inH 2 , 



where n = [(E k — e k )/(X J <\>\ dx)] l l 2 . Moreover, we have dE k Q k ,E k = 0(n~ 2 ) Q k: E k + 
0(n) = 0(n~ l ), and d E Q k ,E = 0(n~ 3 ). For ground states we will drop the subscript 
and write Qe and Re = OeQe- We have Re = Cu~ 2 Qe + 0{n). If we define 
Ci = (Qe, Re)' 1 , we have C\ = 0(1) and Xc 1 > 0. 

The following lemma summarizes the renormalization results near ground states. 
They are proven in |19|] using implicit function theorem. The L 2 -subspace M = Me 



will be defined in Lemma 

Lemma 2.2 Let Y\ = Y, defined in ( |1.10| ), or Y\ = L 2 oc , defined in ( |1.8| ). There 
are small constants uq > and Eq > such that the following hold. Suppose ip G Y\ 
is close to a nonlinear ground state Qe^ iQ with \\ip\\ Yl = n , \\4> ~~ Qe e%B || y < Tn ; 
< n < n , < r < e . 

(1) There are unique small a,9Gl and h G Me such that 

^=[Q E + aR E + h]e i ^ +e \ (2.2) 

Moreover, \\Qe\\y 1 ~ n > a = 0(rn 2 ), h = 0(rn) and 9 = 0{r). 

(2) (best approximation) There is a unique E* near E such that the component 
along the Re, direction as defined by (|2.2| ) vanishes, i.e., there are unique small G K 
and G Me, such that 

^=[Qe, + K] e ^ +9 *\ 
Moreover, E — E* = 0(rn 2 ), = 0(rn) and = 0{r). 
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(3) Suppose E' — E + 7 with | 7 | < r 2 n 2 . By part (1) we can rewrite ip uniquely 
with respect to E' as 

^ = [Q E ,+a'R E , + ti] e^ @+e '\ 
where h! G M E >; h! , a' and 8' are small. We have the estimates 

E + a-E' -a' = 0(r 7 ), h - h! = 0(^ 1 r 7 ), 6 - & = 0(r^ 2 r 7 ). (2.3) 
Notice that n _2 r 7 < Ceq is small. 

2.2 Linear analysis 



We first recall some local decay estimates for e ttH °. The decay estimate (2.4) is 
proved in |||, ^4| using estimates in [§, . The estimate (p.5| ) is taken from |17|, [19 



The estimate (|2.5| ) holds only if we take r — > 0+, not r — > 0— . 



Lemma 2.3 (decay estimates for e Suppose that Hq = —A + V satisfies the 

Assumptions A0-A2. For q G [2, 00] and q' = q/(q — 1), 



e~ itH ° p c <p\\ Lq <c\tr 3 (^) uiy. 



(2.4) 



For sufficiently large r 1; for all k, I, m G {0, . . . , N}, we have 

„-itHo 



lim 



(x) 



-n 



(if + e fc - ei - e r , 



Pc(x) 



-ri 



<C{t) 



-3/2 



L 2 



(2.5) 



L 2 



We now consider the linearized operators around nonlinear ground states. Let 
Q = Qe be a nonlinear ground state with ||Qe|| L 2 small. If we consider solutions 
ip(t,x) of ( |1.1| ) of the form 

4>(t,x) = [Q E {x) + h{t,x)\ e~ iEt , 

with h(t, x) small in a suitable sense, then h(t, x) satisfies 

dth = Ch + nonlinear terms, 

where the linearized operator £ = C E is defined by 

Ch= -i{(-A + V - E + 2\Q 2 )h + \Q 2 h) . (2.6) 

The properties of £ are best understood in the complexification of L 2 (M 3 , C). 
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Definition 2.1 Identify C with M 2 and L 2 = L 2 (R 3 ,C) with L 2 (R 3 ,R 2 ). Denote by 
CL 2 = L 2 (R 3 ,C 2 ) the complexification of L 2 (R 3 , R 2 ). CL 2 consists of 2- dimensional 
vectors whose components are in L 2 . We have the natural embedding 



y.feL 2 



Re/' 
Imf 



eCL 2 



We equip CL 2 with the natural inner product: For f,g G CL 2 , f = g = [%\, 

we define 



if, 9)) 



f-gd 3 x 



{fi9i + hg2) d 3 x. 



(2.7) 



L :;:i 



Denote by RE the operator first taking the real part of functions in CL 2 and then 
pulling back to L 2 : 

r/i 



RE : CL 2 
We have RE oj = id L 2. 



RE 
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(Ref)+i(Reg). 



The operator £ can be naturally extended to an operator acting on CL 2 with the 
following matrix form: 

L_ = -A + V - E + XQ 2 
L + = -A + V - E + 3AQ 2 . 



L_ 
-L, 



where 



(2.8) 



Recall the Pauli matrices 

They are self-adjoint. We have RE £ = £ RE and 

o x £ = £Vi, a 3 £ = -£a 3 , 



"0 


1" 




"0 -i 




"1 " 


1 





, 0- 2 = 


i 


, ^3 = 


-1 



(2.9) 



— -L-i 
L- 



where £* has the matrix form 

We summarize the properties of £ in the following lemma, whose proof is the same 
as that in 0] and Theorem 2.1]. For convenience of notation, we identify L 2 as 
a subspace of CL 2 and make no difference between ip e L 2 and j(^) G CL 2 . 

Lemma 2.4 (spectral properties) Suppose the assumptions AO- A 1 hold. Let Q = 
Qe be a nonlinear ground state, ||<3_e|| L 2 = n, < n < n . Let £ = £ E be defined as 
in (pp. 

(1) The eigenvalues of £ are and ±iu)k, k — 1, . . . ,N, where uj^ = e^ — eo + 0(n 2 ) 
are real and positive. All eigenvalues are simple except which has multiplicity two. 
The continuous spectrum of £ is 



S c = {si: s G R, \s\ > \E\} . 



(2.10) 
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There is no embedded eigenvalue. The bottoms of the continuous spectrum, ±iE, are 
not eigenvalue nor resonance. 

(2) The O-eigenspace is spanned by [q] and [q]. Note £ [q] =0 and [q] is a 
generalized ^-eigenvector with £[$] = — [q\- We denote 

CL 2 , (or S = span {iQ, R}). 

R 

(3) For each eigenvalue iuok, k — 1, . . . , N, there is an eigenvector of the form 
$fc = [-^ fe ], where Uk and Vk are real-valued L 2 -functions satisfying 

L + u k = uj k v k , L_v k = uj k u k , (u k , v k ) = 1. 

Moreover, they are perturbations of (f> k : Uk,Vk = 4> k + 0{n 2 ). §k = is an 

eigenvector of —iuok ■ We denote the combined eigenspaces of ±u>k as 

CE fe = span $ fc } C CL 2 , E k = span 

C R 

(4) The continuous spectrum subspace, H c (£), is equal to 

H c (£) = G L 2 : ((^V, /)) = 0, V/ G S © Ej © ■ ■ • © E^} . 

(5) The space L 2 (IR 3 ,C) ; as a real vector space, can be decomposed as the direct 
sum of N + 2 C-invariant subspaces: 

L 2 (R 3 , C) = S © Ex © • ■ ■ © Eat © H c (£). (2.11) 

For any f and g belonging to two different invariant subspaces, we have the orthogo- 
nality relation 

((aj, g)) = 0. (2.12) 

(6) For any function ( k G E&, k — 1, . . . ,N, there is a unique a k G C so that 

( k = RE a k ^ k - 

Since £ RE = RE C, we have 

C( k = RE ioj k a k $ k , e tc (k = RE e u " h a k $ k . 

(7) By the orthogonality relation ( f2.12j ), any ip E L 2 can be decomposed with 
respect to ( |2.11| ) as 

N 

ip = aR + biQ + ^ RE at k $ k + rj, (2.13) 

k=l 










R 


span < 

R I 


Q 


5 






Uk 









CL 2 
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with 7] e H c (£), a k = c k + id k , 

a = (Q,i?) _1 ((5,Re^), c k = (u k , r fc ) _1 (r fc , Re-0), 
b = (Q,R)-\R, Imip), d k = (u k , v^ 1 ^, haip). 

Note (Q, R)' 1 = a = 0(1) and (u k , v k )- 1 = 1. 

(8) LetM = Eie-'-eEjveHc^). We have L 2 (R 3 , C) = S®M and M = 

For m = 0, 1,2, i/iere zs a constant C > 1 sitc/i £/ia£ ; /or all <p £ M C\ H 2 and all 
itGl, we /iave 



(2.14) 



cr 1 



< e 



< c 



(2.15) 



As in [ID], in order to prove various estimates and make explicit computations, 
we will introduce an L 2 -subspace X and two operators A : X — > X and U : M — > X 
so that 



C 



M 



U-\-i)AU. 



Explicitly, let X be the L 2 -subspace orthogonal to Q: 



X = II(L 2 



2 /TD>3\ 



x 



(2.16) 



(2.17) 



where II is the orthogonal projection which eliminates Q-direction: Uh = h — |qq) Q- 
Let Pm be the projection (not orthogonal) from L 2 onto M according to the decom- 



Pi 
P 2 



position L (K. ) = 5 © M. Pm has the matrix form 
P 1 and P 2 ar e given by (ci = (Q, -R) -1 ) 

P 1 :L 2 -^Q ± , P 1 = id-c 1 \R)(Q\, 
P 2 :L 2 ^R ± , P 2 = id - Cl \Q) (R\. 

Clearly P\R = and P 2 Q = 0. One can check easily that the maps 



where the projections 



(2.18) 



R 



MX 



I 

o n 



: M — > X 



R 



XM 



I 

P 2 



: X 



M, 



are inverse to each other. We now define H = L_ and 

A = [{H 2 + H^mXQ^H 1 ' 2 )] 112 = [H^L+H 1 ' 2 ] 1 ' 2 



(2.19) 



(2.20) 



A is a self-adjoint operator acting in L 2 (IR 3 ), with Q as a 0-eigenvector. We shall 
often view A as an operator restricted to its invariant subspace X. Define 

r A i/2 H -i/2 



U 









4-1/2^1/2 



(2.21) 
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and let 

U = U R MX : M — ► X, ET 1 = RxmUq 1 : X — > M. (2.22) 

Notice that H~ l l 2 is denned only in Q- 1 . We summarize the properties of A and £7 
in the following lemma. 

Lemma 2.5 (Similarity relation) (1) LetlL, A, U and U~ l be defined as in ( [2.171 ), 
( CT ) and (P^2]) ; respectively. Then (g7Tg ) /io/ds in M. // we denote 6t/ P fe A ; 
k — 1, . . . , N, and the orthogonal projections onto the eigenspaces and continuous 
spectrum subspace of A, we have 

UPf = P£UP M , U P c £ = P C A fLP A/ . (2.23) 

(2) The operators U : M ^ X. and U' 1 : X -»• M are bounded in W k,p and L 2 
norms for k — 0, 1, 2, 1 < p < oo, and \r\ < r%. Here r\ > is determined by ( [2.30|) 
later. The operator U — 1 is local and bounded by n 2 , and hence so is [U, i] = [U — 1, i], 
in the sense that 

||[f/,#|| LlnL5 /4<Cn 2 ||0|| L5 . (2.24) 

We have 

U = U + + U_ conj , ET 1 = U* + -U*_ conj , (2.25) 

where conj is the conjugation operator with the Pauli matrix a 3 as its matrix form, 
and the duals of U + and U- are respect to the L 2 inner product (|2.1|) . The operators 
U + and U- are given by 

U± = ^{IiA 1 l 2 H- l ' 2 P l ± UA-^H 1 / 2 ]!), (2.26) 

(U±)* = ^(P 2 H- 1/2 A 1/2 Tl ± IiH l/2 A~ 1/2 Ii). (2.27) 

In particular, U + = 1 + 0(n 2 ), U- = 0(n 2 ). They are not self-adjoint but they 
commute with i and conj . 

(3) For m — 0, 1, 2, there is a constant C > 1 so that 

C- 1 U\\ Hm < \\e- UA 4 Hm < C U\\ Hm , (2.28) 
for all (f) G X n H m and all t G JR. For q E [2, oo] and q' = q/{q - I), 

\\e~ itA P c A n0|| Lq < C\t\<^) U\\ Lq , . (2.29) 
For sufficiently large r% > and for all k, I G {1, . . . , N}, we have 

<C(r 3/2 M L2 , (2-30) 



(*> ~ ri e~ ltA (A -. 1 r P C A n (x) ^ <P 



L 2 ' 

L 2 
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where Oi means ri with lim r _ >0 + outside of the norm. Finally, 
1 



Wl , im - n ; — p c A n0 o 



A — Oi — 2uj k 

1 

Ho-E-Oi- 2uo k 



Q (j>\ , Im - - P C O 02 ) + (n 2 ) > 0. (2.31) 



Estimate ( 2.29|) for A = —A + V was proven in 24| using estimates from 



11 1 . Estimate ( p.30|) for A = (—A + V + m 2 ) 1 / 2 was proven in |17|. Lemma |2.5| is 



a summary of [19|, Lemmas 2.5-2.9]. We omit the proof. 

3 Stabilization regime 

In this section we study the dynamics of the solution when it is close to nonlinear 
ground states. We want to show that the solution ip(t) converges to some nonlinear 
ground state locally as the time tends to infinity. In this time regime, the natural 
decomposition of ip(t) is 

m = [Q E + a{t)R E + C(t) + V (t)} e- iEt+ie W (3.1) 

with respect to a fixed E. Here a(t),6(t) G E, ((t) = (i(t) H hCiv(t), Cfc e E k (C E ), 

k = 1, . . . , N, and rj(t) e H c (£ s ); see Lemma |2^ . Define 



7o" = max lim Im ( o fc 0z, r P c o 0fc<^J, (3.2) 

i<fc,z<Af r^o+ y —A + 1/ + e — efc — e\ — s — ri J 

\s\<sq 

and (recall C\ = (Q, i?) -1 and 7 is defined in ( |1.5|) ) 

D = 6iV| Cl | 7o + / 7o = 0(l). (3.3) 
We will prove the following theorem. 

Theorem 3.1 Assume the assumptions A0-A3. There are small constant uo,Bq > 
such that the following holds. Suppose that the initial data ipo with \\iPoWh 1 <^ 1 *s 
close to a nonlinear ground state QE e %e ° in Lf oc -norm with \\Qe \\l 2 = n — n o> an< ^ 
that in the decomposition ( |3.1| ) of ip with E = E one has 

1 /2 

(j2k=i\\Ck,E \\ 2 L 2^ < |p , \a Eo \<D P 2 , Po<£ n. (3.4) 

Suppose, furthermore, for all E close to Eq with \E — Eq\ < 3-Dpg, the dispersive part 
t]e(0) in the decomposition (|3.1| ) satisfies 

\\e sC VE(0)\\ L5 <n i/5 p(s) 8 / 5 , 
\\e sC VE (0)\\ L2 <A(s) = (l + s)- 1 / 2 p 2 (s) {S " >] 



J 1 o c 
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for all s > 0, where 

P {s) = [pf + N~ l l0 n 2 sy l/2 . (3.6) 

Then there is a frequency with {E^ —Eq\ < 3DpQ and a function 0(i) = — E 00 t + 
0(\og(t)) for t G [0, oo) such that, for some constant C3 > 1 independent of n, 

\\m-QE^\\ L2 <C 3 p(t), (f>0). (3.7) 

loc 

1 /2 

Suppose, furthermore, that (^J2^ =1 ||Cfe,Bollz,2 ) — \Po- Then we also have a lower 
bound 

C^p{t) < \\m - Q Eao e^\\ L2 , (t > 0). (3.8) 

loc 

We call Eq. (|3.5|) the out-going estimates of rj(0). The theorem holds true if they 
are replaced by the following stronger but simpler assumption that 

^0 e Y = H 1 n L 1 ^ 3 ), 11^(0)1^ < pi (3.9) 

since ( p.9| ) implies ( ^.5[ ) by Lemma p75| . Eq. ( |3. 9] ) means that the data ip is localized. 
In contrast, Eq. (pj.5[) only requires the data to be "out-going" in some sense. This 



will be useful when we prove Theorem using Theorem |3j] in section 4 



Our strategy of proof is as follows. For each T > 0, we choose Qe{t) to be the 
best approximation of ip{T) given by Lemma |2.2| . We will prove estimates for the 
components of ifj(t) in the decomposition (|3.1|) with respect to E — E(T) for t 6 [0, T]. 
We will then use a continuity argument to show that E(T) can always be chosen and 
we have uniform estimates as T — ► 00. It then follows that E(T) converges to some 
E^ close to E(0) as T — >• 00. 

3.1 Equations 

Let Q = Qe be a fixed nonlinear ground state with frequency E near eo, and i? = 
Re = 9eQe- We write the solution ip{t,x) of (|1.1| ) in the form 

^(t, s) = [Q E (x) + a(t)R E (x) + h(t, x)) e - iEt+ie{t \ (3.10) 

where a(t), 6{t) G K and h(t, .) G M E . M s is defined in Lemma ^]4| (8). Substituting 
the ansatz (|3.10|) into ( |1.1|) and using CiQ = and CR = —iQ, we get 

d t h = Ch + i~ l (F + 6(Q + aR + h)) - aiQ - aR. (3.11) 

Here 

F = XQ(2\h a \ 2 + h 2 a ) + X\K\ 2 h a , h a = aR + h. 
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We want to choose a(t) and 9{t) so that h(t) G M, that is, /i(0) G M and 



L (F + 6{Q + aR + h)) - aiQ — aR G M. Since M 



a 1 



a(t) and satisfy 



(o, lm(F + Oh) - dR^j = 0, 
R, Re(F + 6(Q + aR + h)) - aQ) = 0. 



Denote c\ = (Q,R) 1 . We have 



a = ( Cl Q,lm(F + 9h)), 

= -[a + ( Cl R, Re F)] ■ [1 + a( Cl R, R) + (^R, Re h)]' 1 



(3.12) 



Eq. ( |3.11|) for h becomes 

d t h = Ch + P M F all , F all = r\F + 9{aR + h)). 

We decompose h(t) with respect to the spectral decomposition ( 2.11]) , 

h(t)=((t)+T)(t), C = Ci H + Cjv, 

where (k G E&(£), = 1, ... ,N, and r\ G H c (£). For each Ck{t) we associate a 
function z k (t) G C by writing ( k = ~Re(z k ) u k + lm(z k ) iv k - In other words, Cfc(^) = 
RE z k (t)<& k . H we define = \{u k ± v k ), we can write 



Cfc = Re(zfc) u k + Im(^ fe ) it> fc = z k u\ + ^Uj. 
Projecting (|3.12|) to E fc (£) we get 



(3.13) 



Zk 



-iu k z k + Re Fan) +i(uk, ImF ( 



all ) 



(u k ,Faii). 



From the linear part, we identify the phase factor of z k as — uo k . Hence we define 
Pk{t) = e luJkt z k (t), which has the same magnitude as z k but with no strong oscillation. 
p k (t) satisfy 



e^ t p k (t) = (ui,F all )-(u k -,F, 



all) 



1 { («+ F) + (Ufc, F) + [(«+ /») + ( % -, &) + K, i?)a] . 



Also, projecting (|3.12|) to H c (£) we get d t r} = £r]+ P c c F aU . 
Summarizing, for 

^(t) = (Q + a(t)R + h(t)) e - iEt+ieit \ 

h = ( + r), C = CiH 1- Cat, Cfc = ZkU k + z k u k , z k 



Pk, 
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we have 

a=( Cl Q,lm(F + 9h)), 

ie^h = (u+,F) + (u^,F) + [(u+,h) + (u- k ,h) + (u k , R)a] 9, (3.14) 
d t r) = Cr]+ P c £ r x (F + 9(aR+h)), 

where c\ = 

F = \Q(2\h a \ 2 + h 2 a ) + \\h a \ 2 h a , h a = aR + h, (3.15) 

e = F e = - [a + ( Cl R, ReF)} ■ [1 + a( Cl R, R) + (aR, Re h)]' 1 . (3.16) 



This is a system of equations involving a, z k and rj only. Note that 9 enters ( p.14 ) 
only via 9 = Fg. It will appear in the form e td when we integrate 77. Hence we do not 
need estimates of 9 for the proof. 

For convenience, we will use the following convention. 

Convention 3.1 For k — 1, 2, • • • , N , denote 

u-k = -w k , z_ k = z k , p_ k = p k . (3.17) 
We have z k (t) = e~ Wkt p k (t) for both k > and k < 0. We also denote 

n = {±1,... ,±N} . (3.18) 

3.2 Decompositions of F, a and 77 

Most quantities in our system of equations are strongly oscillatory. It is necessary to 
identify their oscillatory parts before we can estimate. In this subsection we identify 
the leading oscillatory terms of a and 77, and decompose F according to order. We 
will treat z k and a again in §3.3. Note that 

Q = 0(n), R = 0(n- 1 ), c x = 0(1), M + = ^ + 0(n 2 ), u ~ = 0(n 2 ). (3.19) 

In fact, since ||(5e || L 2 = n and \E — E \ < 3Deln 2 , we have ||(5b|| L 2 = [1 + 0{e^)\n 
and Q — (1 + O(el))n<f) + 0(n 3 ). We will prove that 

\z k (t)\ < Cr 1/2 , \a{t)\ + \\r](t)\\ L 2 <Cr\ ast^oo. (3.20) 

loc 

Hence the main term in h a = aR + ( + 77 is (. Therefore, the main part of F, defined 
in ( |3.15| ), is 



F 1 = AQ(2|C| 2 + C 2 )- (3-21) 
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3.2.1 Decomposition of a 



We now identify the main oscillatory terms of a(t). Recall from (|3.14j ) that a 



(ciQ,lmF + Oh), c\ = We shall impose the boundary condition of a at 

t = T, which is in fact a condition imposed on the choice of E(T). Hence we use the 
following equivalent integral equation: 

a(t) = a{T) + j (ciQ,ImF + Oh){s) ds. 

The main term of Im(F + Oh) is ImFi = lmXQ( 2 . Thus the main oscillatory terms 
of a(t) are from the integral ft ds with 

N 

A< 2 > = ( Cl Q,lmXQ( 2 ) = ( Cl XQ 2 ,lm ^ (kCl)- 



k,l=l 



Since 



ImC fc C; = lm(z k u k + z k u k )(ziuj + Ziu x ) 

= \m(z kZl )(uiui - u^ul) + ImizkZ^iu+uJ- - u k u+), 



we have 



N 



£(2) = { Qkl l i m ( ZkZl ) + aM 2 lm(z k zi)} , 

k,l=l 

where au,x = (ciXQ 2 ,(u k ul - u k uf)) and a k i :2 = {ciXQ 2 ,u k uJ - u k u~l) are real 
constants bounded by n 2 . We can integrate by parts A^ to get: 

/ A (2) ds = Im \~] / a k i A (z k z t ) + a Mt2 (z k zi) ds 
Jt k,i=i jT 

= Im V / e-^^aM^pm) + e^-^ s a kl ^p k p{) ds 
k,i=i jT 

= Im ^ [ ia ki,3Z k zi + ia kiA z k z^\ T - I A 2jTmd (s) ds, (3.22) 

kl=l ^ T 



where 

0-kl,3 = (^k + Ml) 1 a k l,l, <Xkl,A = 5 k {iO k — LUi) 1 a k i y2 , (3.23) 

N f d d ^\ 

A 2 , rmd = Im J2\ e- l{ulk+u)l)s ta kl , 3 -(p kPl ) + e^-^'ia^—ipm) \ . (3.24) 

k,l=l ^ ' 
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Here we put 5 l k in the definition of a^ i4 to impose a k k,4 = 0. We have the factor 5 l k 
in dki,4 since 1m z k zi = if k = I. Also note that 0^,3 and au,4 are real constants 
bounded by n 2 . The first part in ( |3.22| ) can be rewritten as aP> (t) -a<$(T), where 



N 



a (2) = Im ) j ia k i^z k zi + ia HA z k zi = ^ a M z k zi. (3.25) 
k,i=i k,ien 

Here we have used Convention |3.1| . The constants a k i = \a k i^ if k and I have the 
same sign; a k \ = \a k i^ otherwise. In particular, a k \ are real constants bounded by n 2 . 
aS 2 \t) contains the main oscillatory part of a. We denote the rest of a{t) by b(t), 

a(t) = a {2 \t) + b(t). (3.26) 

Thus b(t) = a(T) - a< 2 >(T) + J^b(s) ds with 

6 = (ciQ, lm[F -F 1 + 6h))- A 2)Tmd . (3.27) 

As we will see later that b(t) is smaller than d^ 2 '(t). However, b(t) is the main part 
of a(t): We have b(t) < p 2 (t) while a^(t) < n 2 p 2 (t). 

3.2.2 Decompositions of F 

Recall ( ETT5D , 

F = XQ(2\h a \ 2 + h 2 a ) + X\h a \ 2 h a , h a = aR + h. 



In view of ( |3.20| ), we decompose h a = ( + bR + (rj + a^R) and decompose F as 

F = F x + F 2 + F 3 + F 4 + F 5 , (3.28) 

where 

F 1 = Ag(2|C| 2 + C 2 ), 

F 2 = 2\QRb(2( + C) + 3XQR 2 b 2 + \(( + bR) 2 (( + bR), 
F 3 = 2\QRa®(2( + (), 

F 4 = 2\Q[(( + C)r/ 2 ) + Cn { % (3.29) 
F 5 = 2XQ[(C + C)V {3) +CV {3) } 

+ \Q [2\r] a \ 2 + r/ 2 ] + 2\QRb(2 Va + %) (rj a = rj + a®R) 

+ X{aR + h) 2 {aR + h) - A(C + bR) 2 {( + bR). 

Here F\ consists of terms of order nz 2 ; F 2 , F3 and F4 consist of terms no smaller than 
n 2 z 3 ; and F5 higher order terms. 
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3.2.3 Decomposition of rj 

We now identify the main terms in rj. We first recall from (|3.14j) that 

d t rj = £rj+ Pfr^F + 6(aR+( +77)]. 

Using Lemma p] that C = TJ- x {-i)AU on H c (£) and UP C C = P A U, we have 

d t Urj = -lAUr] + P A Ui'^F + 6(aR + ( + rj)] 

= -iAUr} - idUr] + P C A Ui'^F + 6(aR + £)] - P C A [[/, i]^. 

Here we have used the commutator [U, i] to interchange U and i so as to produce the 
term iOUrj. This term is a global linear term in Urj and cannot be treated as error 
(however [U, i}U~ l 6r] is a local error term). We can eliminate it by introducing 

fj(t) = e im U7](t), 6(t)= [ F e (s)ds. (3.30) 

Jo 

We have fj(0) = Ui](0) and 

d t fj = -iAfj + e ie P C A Ui' 1 [F + 6(aR + ()] - e ie P C A [U, i]9r}. 
Hence fj(t) satisfies the integral equation (using e~ lAt rj(0) = Ue tc r](0)) 

fj(t) = Ue tc r]{0) + f e - iA{t ~ s) P A F v {s) ds, (3.31) 
Jo 

F v = e iB U%- x [F + 6(aR + C)] - e ie [U,t]6r]. (3.32) 

Since U and U^ 1 are bounded in Sobolev spaces by Lemma |27|, and 

v (t) = U^e^fKt), (3.33) 

for the purpose of estimation we can treat rj and fj as the same. 
To identify the main term of fj, we decompose F v as follows, 

F v = F Vt2 + F n $, 

F v , 2 = e ie m~ 1 F 1 , (3.34) 
F v , 3 = e %e Ui~ l [{F - F x ) + 6{aR + C)] - e* e [U, i]6rj. 

The leading part of fj is from F Vt2 - Recall F x = \Q(( 2 + 2((), and U — U + + conj \J_ 
with U + and U- commuting with i and conj , see ( |2.25|) . Hence 

F v ,2 = e^Ut-'XQiC 2 + 2CC) 

= e i6 i-\U + - conj U_)\Q(( 2 + 2(C) 

= e ie i~ l {U + \Q(( 2 + 2CC) - U^XQ(C 2 + 2(0} • 
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Substituting ( = Y^aLi (i and d = ziu\ + z i u l , we have 

N 

F v , 2 = eV 1 {u+XQ(( k (i + CkCi + GO) - U-\Q((k(i + GO + GO)} 
k,i=i 

In the last line, Convention p. 1| is used. In particular, z\. = z\u if k < 0. The functions 
$fcz are defined as follows. For k < 0, denote = and -u^ = ujL. We define 

k,l>0: $ ki = U + \Q(u~lu+ + u£uf + % u~f) 

- U_\Q(u- U ; + u+itf + v^uf), 
k > 0, 1 < : $ fci = U+XQiu+ul + + UfcU,") 

- f/_AQ(w fc w+ + « + u^uf ), (3.35) 
£;<0,Z>0: $ w = $; fe , 

fc, Z < : $ M = U + \Q(il^uJ + m£ + m+) 

- U_\Q(ulut + U+lif + %M+). 

Since = 0{n 2 ) and f7_ = 0(n 2 ), we have $ H = 0(n 3 ) if A;, I < 0, $ w = 0(n) 
otherwise. 

We now integrate F Vj2 - Using -Zfc(s) = e~ lu}kS pk(s) (see Convention |3.1| ) and inte- 
grating by parts we get 

f e~ iA ^ P A F v>2 ds = T f e-^^e^z.zt P A *« 
Jo Jo 



-E 



as 



where the summation is over k,l £ Q and 
%z = lim 



(Jfe,Z e fi). 



(3.36) 



■r^o+ A — u;^ — u i — ri 
Here we add — ri since A — lo^ — lo\ may not be invertible. We take r — > 0+, not 
r — ► 0—, to get the decay of the second and the last terms above, see Lemma |2]5| (3). 
Recall cjfc ~ efc — eo > for k > and to^ = —uj\k\ if A; < 0. Since A has a spectral 
gap ~ |eo|, A — ujk — toi is not invertible only if k > and Z > 0, by Assumption A2. 
Hence rj^i e L 2 unless both and Z are positive. 
We denote the main term above by fj^ 2 \ 



E 

k,ien 



(3.37) 
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and denote the rest of 77 by r]( 3 \ We decompose fj^ = r/^ + • • • + rif 1 where 
fjf* = Ue tc r](0), 

# = " E e- iAt (e i9 z kZl )(0)7} kl = -e~ iAt fP (0), (3.38) 

^ 3) = - E f^^-" ^^w*))^, 

Men 70 dS 

# = f e-^(*- s ) P C A F, i3 da. 

We have 

Correspondingly, we define 

^ = u-^fP = u- 1 E 

Men (3.39) 
^(3) = C /-i e -i^(3) > „(3) = U-le-Hfjf), (j = 1, . . . , 4). 



Lemma 3.2 Suppose, for a fixed time, for some p Cn < n 0; 

IIQ||=«, IMI^nL* « 1> NIl? <«> 

loc 

max \z k \ < p, \a\ < Cp 2 . 



Denote 



X = p 2 \\v\\ Lf +n\\v\\% +|h 3 || Ll 



J loc ^loC 11 ' " "^loc 

|2 , ||„3| 



X = np\\ri\\ Lf +n\\ V \\ L 2 + \\v 6 \\ L i 

loc loc " " -^loc 



(3.40) 



(3.41) 



We have 

\\F4lI < np A + np\\ V W\\ L2 +X, 

loc ^l oc 

\\F 3 + F A + F 5 \\ L} <n 2 p 3 + X, 

loc 

<P 3 + X, \\F\U <np 2 + X, 

loc loc 

\\F e \\ Llc <p 2 + n- l X, (3.42) 
m ax \p k \ < np 2 + X, \b\ < np 3 + nX, 

k 

\\ F nA\ifi/W ~ P 3 + n P IMIl? oc + n W r )\\ll c nL5 + lh 3 Ls/4 ni i +"V WvWl* > 

II^IIl5/4 < V + ^pII^IIl? + w IMIl? nL 5 + ||^ 3 || L 5/4 + «V hlLs • 
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Proof: By fl3.4U| ) we have |a (2) | < n 2 p 2 , \b\ < p 2 and ||^ (2) || L 2 < np 2 . Also recall 



( |3.19| ). Hence for F 5 defined in ( |3.29| ) we have 



11*5 |Li <np\\v (3) \\ L 2 +n(\\r)\\ L i + up 2 )* + p 2 (\\r)\\ L? + np 2 ) + \\ V 3 \\ 

loc 11 11 -^loc \ loc / \ loc / -^loc 

<np 4 + np\\r]^\\ L2 +p 2 \\v\\ L 2 +n\\ri\\ 2 Lf + \\v 3 \\ L i 
= np 4 + np\\rj {3) \\ L2 + X, 

loc 

where X is defined in (|3.41| ). Combining 7p 2 ) + 77 ( 3 ) = i] in the definition (|3.29|) of F<± 
and F 5 , we can use the same argument for ||F 5 |Li to show 

loc 

ll*4 + *5|Li <np A + np\\ V \\ L? +X<np A + X. 



By their explicit form in (|3.29|) , we have H-F3L1 < n 2 p 3 , ||F 2 |Li ~ \bz\ + b 2 /n + 

loc loc ' 



(z + b/n) 3 < p 3 , and ||-Fi|| L i < np 2 . Hence 

loc 

\\F 3 + F 4 + F 5 \\ L} <n 2 p 3 + X, \\F-F 1 \\ L} < p 3 + X, \\F\\ L} <np 2 + X. 



From definition ( |3.16| ) of Fq we have 

\H < |a| + n- 1 \\F\\ L} <p 2 + n\np 2 + X) < p 2 + n~ l X. 

loc 

By we have 

\Pk\<\\F\\ Ll +(p+\\v\\ L? )\Fe\ 

loc loc 

<np 2 + X + (p+ |MU )(p 2 + n l X) < np 2 + X. 

loc 

By §M, 



\b\ < Cn \\F - Fx\\ L i +\( Cl Q,lmh)\-\F e \ + \A 2irmd \. 

loc 

Note that {(ctQJmh)] < n 3 \\h\\ L i since Imh _L R and Q = Cn 2 R + 0{n 3 ) by 

loc 

Lemma |2.1| . Also |A 2)rm ^| < n 2 pmax^ \p^\ by definition (|3.24|) of A<i )Tm d- Thus 
\b\ < n{p 3 + X) + n 3 (p + |M| L2 )(p 2 + n- 1 ^) + n 2 p{np 2 + X) < np 3 + nX. 

loc 

Let r = 1 or r = 5/4. The estimate for ||F — -Fi|| ir is the same as that for 
\\F — Fi\\ L i except for the non-local term Xrj 2 fj and ||w/ 2 ||rr < n \\r)\\ \% nL5 . Thus 

loc loc 

\\ L r <p 3 + np \\v\\ L? + n \\ V \\ 2 L? nL5 + \\ V 3 \ 



loc 



From the definition ( |3.34j ) of F Vj 3, the boundedness of U, and the estimate ( |2.24p for 



[U, i] , we have 

ll^,3|| L r < II* " F 1 \\ Lr + {FeKn-^al + max \z k \) + n 2 \F e \ \\ V \\ L5 
<p 3 + np\\r]\\ L? +n\\r]\\l ? nL5 + \\r] 3 \\ 



+ (p 2 + n- 1 X)(p + n 2 n , inL5/ 
< p 3 + np\\r]\\ L 2 +n\\r]\\ 2 L 2 nL5 + \\v' 3 \\ L r + n 2 p 2 U nn, -, . 
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Here we also used ||^ 3 || L i < ll 7 ? 3 !!^- Since \\F V \\ ,5/4 < np 2 + ||i^,3 || L 5/4, the estimate 

loc 

for ||Fj L5/4 in (|3~22j ) follows. Q.E.D. 

3.3 Normal forms for equations of bound states 

Recall that z k (t) = e~ lulkt p k (t) and a(t) = a^ 2 '(t) + b(t). Since many terms in the 
equations of p k and b are oscillatory and only contribute to the phases, we now de- 
rive the normal forms for the equations of p k and b, where those terms are removed. 



Throughout this subsection, Convention |0| for uj k , z k and p k with k < is under- 
stood. 

Lemma 3.3 (Normal Form) Suppose, for some p(t) <Cn< uq, 

\\Q\\=n, IMi)|| L2nL 5«l, \\v(t)\\i^<n, 
max |* fc (t)| <p(*), \a(t)\ < Cp 2 (t). (3 ' 43) 

fe=l,... ,N 

There are perturbations q k (t) of p k (t), k = 1,... ,N, and (3{t) ofb(t), to be defined 
in ( |3.96| ) and ( |3.103| ), satisfying 



such that 



q k -Pk\<Cnp 2 : \b- P\<Cnp A + n z p\\r ] \\ L2 , (3.44) 



q k = ^ D kMi\ 2 Qk + iJkqk + 9k, (3.45) 

1=1,... ,7V 



P= B kl \z k \ 2 \ Zl \ 2 + g b , (3.46) 

l<k<l<N 

where J k (t) are real functions bounded by Cp 2 to be defined in ( [3.931 ). The constants 
D k i are complex and bounded by n 2 , with their real parts given by ( |3.72|) satisfying 

-57 + n 2 < ReD kl < Cn 4 , Re D kk < - lo n 2 . (3.47) 

The order one constants 70 and 7^ are defined in ( p..5|) and (|3.2|) . The constants B k i 
are real and B k i = — | c\ ReD k i + 0{n A ). Moreover, 

\9k\ < np 4 + n 4 p 3 + np\\r]^\\ L2 +X, (3.48) 



\9b\ <np 5 + n 4 p 4 + np 2 \\rj ( - 3) \\ r2 + n 3 p 2 \\r]\\ T 2 +n 2 \\r]\\ 2 r2 +n\\r] 3 \\ Tl , (3.49) 
where X = p 2 \\r]\\ L 2 +n\\r]\\ 2 -2 + II^IIl 1 is defined in (|3.41|) . 

loc loc loc 
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PROOF: We first observe a consequence of Assumption A3: Suppose k , k±, . . . , kj G 
{1, . . . ,N} with j < j max = 3 satisfy 

u ko ±uj kl ±---± u kj = 0. (3.50) 

Then j is odd and, after a relabeling, the two sets (with multiplicities) |fc , ■ ■ ■ , kj-i j 

and , ... , are the same. The reason is as follows. Recall cu k = e k — eo + 

0(n 2 ) > 0, and hence the expression in ( [3.50|) is approximately equal to 

(e fco - eo) ± (e kl - e Q ) ± • • • ± (e kj - e ). (3.51) 

For the expression in ( p. 50|) to be zero, one of the uo km must be with a negative sign. 
Hence two eo's in ( |3.51| ) cancel and (|3.51j) is the sum of j e k s minus the sum of 
another j e k s. Since j < j max , Assumption A3 then implies that they can be divided 
to two equal sets. Since there are (j — 1) eo's and all fc m 's are positive, j must be 
odd, these eo's cancel each other, and the other e^'s form two equal sets. 
Part 1. Excited states. Recall from ( [3.14|) the equations for p k : 



Pk 



"V^ [«, F) + (u~, F) + [(u+, h) + (u-, h) + (u k , R)a] F e ] , (3.52) 

where h = d + ■ ■ ■ + (n + r]. The nonlinear terms F and Fg are given in ( |3.15|) -( j3.16|) 
and F = F\ + ■ • ■ + F§ is decomposed in (|3.29|) . We shall derive the normal form of 
( |3.52| ) using integration by parts for those terms no smaller than n 2 p 3 . 

Step 1 Integration of terms of order nz 2 . 

Recall that the main part of F is F\ = AQ(2£C + C 2 )- The main terms on the 
right side of ( 3.52j ) are among the first two groups of terms involving F, 



- 1 e i »« t {(u+,F 1 ) + (u k -,F 1 )}. 



We denote the rest by R k 

Rk 



F l ) + (u^,F-F 1 ) 
+ [«, h) + (u k ,h) + (u k , R)a] F e 



and rewrite (|3.52 ) as 



Pk 



- 1 e^{( M +,F 1 ) + KT,F 1 )}+ J R fc . 



(3.53) 



(3.54) 



(3.55) 



Substituting F 1 = Xw, m =i ^Q( 2 0Cm+0Cm) with Q = Ziuf+zty and using Convention 
3TT| , we get 

rV^{(«+ F 1 ) + (u^,F 1 )} 



N 



E 

l,m=l 

E 

l,m £ Q 



"V^* {(XQu+, 2C,Cm + CtCm)) + (AQ% , 2CiCm + CiCn)} 



Im iuj k t 
u k ° ^l^mi 
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for some purely imaginary constants c l k m bounded by n. The phase factor of a typical 
term e lUkt z\z m is io^ — uji — u m , which is nonzero by the observation for ( 3.5(J| ). Hence 
we can define 



-ic lm 



Pk,i-= T k ~ e^z lZm , (3.56) 

< ' f,h (ill f,l... 



0J k - OJi - u, 
■lm 



P M e= E " e^-^^ipWrn+PiPm), (3.57) 

Lm 6 s2 



and we have 



E ^ e ^ fc< Z l Z m = Pk,l- + Pk,le- 
l,m € f2 

Because Pfe.ie are of order n 2 z 3 , we need to extract their leading terms. Substitute 
( |3.55| ) into (|3.57 ) and collect terms. The leading terms in pu,ie are cubic polynomials 



in Zk and Pk\ e can be rewritten as 



Pk,le = i'^dky z l z m z j+9k,l, (3-58) 

l,m,n G £7 



9k,i 



Y\ — e^-^-^iRiPn+ptRn). (3.59) 

' * l,U Mi I.) 



UJk — UJl — (jJ m 



Here dj^f are real constants bounded by n 2 , and the error terms gk,i are bounded by 



\9k,i\ < C max \c k m \ \pi\ \R m \ < Cnp max \R m \. (3.60) 

k,Lm m 



We conclude 



Pk ~ Pk,i- = Rk+ E i le%UJkt rf M z ^ z ™ z i + 9k,l- (3.61) 

l,m,n G f2 

Step 2 Integration of terms of the form nzrf 1 ^. 

We now integrate i" 1 e* a;fc *('u^, F 4 ), which is part of the first term of R k . Recall 
F 4 = 2AQ[(C + C^ (2) +(V {2) ] with r/ 2 ) = U' 1 Ei, m en z i z ^i m - Denote 

rf = } J ZiZmfjim, (3.62) 
l,m £ S7 



where r]i m are defined in ( |3.36|) with ||?7i m || i 2 < n. Recall Lemma |275| that [/ 1 

loc 



([/+)* — (£/_)* conj and both U+ and £/_ commute with i and conjugation. Then 
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rjW = U- X rj = (U+Yrf - {UJ)*ff and 

(«+ F 4 ) = (2AQu+(C + C), V 2) ) + (2AQu+C, #) 
= (2AQ«+(C + C) ) [(^ + )V-(^-)V]) 

+ (2AQu+C, W + yf? - {U-M) 
= (U + [2\Qut(C + 0) ~ U42XQu+C], rf) (3.63) 

+ (U + l2\Qut(]-U42\Qut (( + ()}, fj'). (3.64) 

Substituting ( = £\ + an d V = Ei.me!) 2 ' 2 *! we can wri te 

r 1 e <w * t (u+F 4 )= ]T 4 mj 'e^*z z z m% , (3.65) 

l,m,j £ f2 

for some coefficients g^" 1 - 7 bounded by n 2 . The phase factor of a typical term e itJfei z\z m Zj 
is co>fc — co>z — uj m — ojj. By the observation for ( |3.50| ), it is nonzero unless one of \l\, 
\m\, \j\ is k and the other two are the same. In this exceptional case e luJkt ziz m Zj is of 
the form e luJkt z k Z[Zi = \pi\ 2 p k . For fixed I > 0, there are six such terms if I ^ k, and 
three terms if I — k. We denote the sum of their coefficients as D k i, 

D kk = d kk ^ k) + d k k { ~ k)k + d[7 k)kk . (3 ' 66) 

The total of these zero-phase-factor terms is YliLi D k i\pi\ 2 p k - The other terms can be 
integrated. Define 

•jlmj 

Pk,2-= E — e^ziZmZj, (3-67) 

E—id lmj d 
k e <K- Wl -^- Wj ) t {p ) (368) 
uj k - ui - u m - Uj dt 



We have 



d 



£ 4 mJ ^ = ^ £^M 2 Pfc + ^(P*,2-) + Sfc,2, ( 3 - 69 ) 



and 



\9k,2\ ^ (maxId^Dp^axIpjl < n 2 p 2 max 1 . (3.70) 

j j 

We now compute ReD k i. We want to collect terms of the form Ce luJkt ziZiz k 
from ( |3.63| )-( p.64 ) with ImC 7^ 0. Hence we only need to consider those fji m with 



Im^ m 7^ 0, i.e., l,m > 0. The only terms from rj' with resonance coefficients 
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are ^2 m j >Q z m Zjfj m j, which are of the form z m Zj, with two bars. Hence the inte- 
gral in ( |3.64| ) does not contain z\ZiZ k and is irrelevant. From the integral in ( 3.68| ). 
we want to choose z\z k from rj', i.e., >0 z m Zjfj m j, and choose zi from ( or ( in 
U + [2XQul(C + C)]-U-[2XQu+C]- (Recall (Q that (f,g) = J fgdx.) For the second 
part we get Zi<&ki where 

$> kl = U + 2XQu+(u+ + Ui) - U_2\Qu+uY. (3.71) 

Since fj mj = fj jm for m,j > 0, the terms with z t z k in ^2 mjj>0 z m Zjfj mj is ( 2 - S l k )z k Zifj kl . 
Therefore, also using the definition (|3.36|) of r} k i, 

ReD kl = Rei~ 1 J $ w (2 - 5 l k )fj k idx 

= -{2 - 5 l k ) Im U kU P C A . (3.72) 

V A - uj k - uji - Oi J 

Recall 

Q = (l + O(e 2 ))n4 >0 + O(n 3 ), u\ = 4> k + 0(n 2 ), u k = 0(n 2 ), 

(3 73) 

U+ = l + 0{n 2 ), U- = 0(n 2 ). 
Hence, by ( |3~35|) and 

$w = AQ0 fe 0/ + 0(n 3 ), $ H = 2AQ0 fc 0, + 0(n 3 ), (3.74) 

and we have 

Re D kl = -2(2 - S l k )(l + o(l))n 2 Im (V^, P c A 0o0fe^ J + 0(n A ). 

(3.75) 

By Assumption A2 and Lemma |2~5| (3), we have — ReD k i < bn 2 ^Q and ReD k i < Cn 4 



for all k, I > 0, and - ReD fcfc > 7 n 2 . We conclude (|3.47| ). 

Step 3 Integration of other terms. 

We now integrate other terms in (|3.61| ) no smaller than n 2 z 3 . We first consider 
the first group of terms in R k defined in fl3.54j ). Using F — F\ = F 2 + F 3 + F 4 + F 5 
and removing i' 1 e luJkt (u k , F4), we get 

{(u+,F 2 + F 3 ) + (u k ,F 2 )}+g k>3 , 

where 

9k,3 = i- l e™ kt {(u+, F s ) + (u k ,F 3 + F 4 + F 5 )} . (3.76) 
We have, using u k = 0{n 2 ) and the explicit form of F 3 + F 4 , 

\9kA^\\ F A\L} +n 2 p ; 3 + F4|| L i <\\F h \\ L} +« 2 (n 2 p 3 ). (3.77) 
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We now consider the second group of terms in R k , see (I 



-V^ [(«+ h) + («" /i) + (u fc , i?)a] 



(3.78) 



The only relevant terms are i l e Wkt [(u k , () + (u k , £)] We move the other part to 
error term, 



g k , 4 = rV^ [{u+, 77) + (^, 77) + (u k , R)a] Fg. 



(3.79) 



Note that (u k , R) = 0(n) since u k -L Q and i? = Cn Q + O(n) by Lemma EO 



(3.80) 



Also, |(Mfc,^) + (u k ,T))\ < Cn 2 \\r]\\ r 2 since Ren _L t^, Im?7 _L and the differences 

loc 

between u k , v k ,u k , u k are bounded by n 2 . Hence 

|<? M |<(n 2 NU +n/)|F,|. 

loc 

Using C = Y.i z i u t + we § et 

feC) + K-,C) = E£i (c\zi + cizi 

where C\ 

0{n 2 ) and = 0(n 2 ), see (|3.73|) , we have 

C l k = 6 l k + 0(n 2 ), C l k = 0(n 2 ). 
Thus, terms no smaller than n 2 p 3 in (|3.78|) 



( u t > + K > u * ) and C fc = K > + K > Because u+ = <p k + 



are 



9,2 1 



where Fg A and Fg^ are leading parts of Fg fl3.16| ), (recall a = 6 + a^ 2 ^) 

F M = -[a+ { Cl R, Re F x + F 2 )] • [1 + { Cl R,R)b}- 1 , 
F 6:2 = -[b+(c 1 R,ReF 1 )], 

and ( 3.78|) is equal to g k> 4 + (P-81 ) + g kj 5 where 



(3.81) 



(3.82) 
(3.83) 



9k,5 = i- 1 e iUkt \z k (Fe-Fe, 1 )+ 1 E l {C{- 5 l k ) Zl + C\z x {F -F 6j2 ) 



(3.84) 



Since \F e - F e<2 \ < \F g - F 6il \ + \Fg tl - F 6;2 \, we have 

\9k,5\ ^ p\Fg- Fg tl \ + n 2 p \Fg A - Fg, 2 \ . 
Recall ( gyp Fg = - [a + ( Cl R, ReF)) ■ [1 + ( Cl R,R)a + (c^Reh)]' 1 . Hence 



Fg - F 0A = - ( Cl R 7 Re F 3 + F A + F b ) ■ [1 + {c x R, i2)6] _1 

_ [a + {ctR, ReF)] ■ [( Cl R,R)a {2) + (dR,Reh)} 
+ [1 + (dR, R)a + (cxR, Re ft)] • [1 + {c^R, R)b] 



(3.85) 
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Since R = Cn 2 Q + 0{n) and Reh J_ Q, we have |(ci_R, Re/i)| < n ||/i||ri • Hence, 



iFe-F^l^n-MlFa + ^ + FslU + |a| + n" 1 \\F\\ L} (p 2 + n N , t || £ i 

loc locj loc 

< n-^nV 3 + X) + [p 2 + n-^np 2 + X)](np + n \\rj\y ) 

loc 

< np 3 + n^X. 



Here we have used Lemma 3.2. We also have 



F 9 ,i - F 0i2 = -[a< 2 > + [ Cl R, ReF 2 )} • [1 + {c x R, R^}' 1 

+ [a+( Cl R, Re F 1 + F 2 )} ■ [1 + ( Cl R, R)^- 1 ■ ( Cl R, R)b. 



(3.86) 



Hence 



\F 9 ,l - F efl \ < (nV + n- 1 \\F 2 \y ) + (p 2 + n" 1 Hi^ + F 2 || £l )n 2 p 2 . 



Since ||-Fi|| L i < np 2 and ||i<2|| £ i < p 3 , we have \Fgi — F e ^\ < n 2 p 2 + n 1 p 3 . Hence 

loc loc 

gkfi is bounded by 

Ifl'Jfe.s | %p\Fe~ Fe t i\ + n 2 p \F e ,i - Fg i2 | 

< p(np 3 + n~ l X) + n 2 p(n 2 p 2 + n~V) 

< np 4 + n 4 p 3 + n- x pX. (3.87) 



Summarizing, we can rewrite eq. ( 3.61|) as 

N 



d 

~dt 



{Pk -Pk,i- ~Pk,2-) = ^2 F>ki\zi\ 2 Pk + Rk + ^2,9k 



i=i 



i=i 



where 



Rh = i 



-V^j 47 W, + (ut,F 2 + F 3 ) + (u- k} F 2 ) 

K l,m,j £ f! 



(3. 



- Zk l + ( Cl R,R)b [a + ^ ReF ^ F ^ 

~ Ei \(C l k - 5 l k ) Zl + C l k z\ [b + (dfl,ReF a )] }■ (3.S9) 



We now integrate i2fc. Denote 

B = ( Cl R,R)b, \B\<Cn- 2 p 2 . 
Note that is of the form 



te 



[n^b 2 + bz + n^bz 2 + z 3 + n~ x z k z 3 \ 



(3.90) 



(3.91) 



31 



in the sense that 

Juj k t 



Rk 



l + B 



Zj 2 



Zj 2 Zj 3 



(3.92) 



+ Yl f( B ) z h z j2 z js+ Yl f( B ^ n lzkZ K' 

where f(B) are polynomials in B with real coefficients bounded by one. We have 
omitted their dependence on the summation indexes. The phase factors of the above 
summands are 

CJfe, UJ k ± U\ h \, U) k ± UJ\ h \ ± LU\ j2 \, U) k ± U\ h \ ± LU\j 2 \ ± LU\ j3 \, ±UJ\ji\ ± W|j a | ± W|j-3|, 

respectively. By the observation for ( |3.5U| ), zero phase factor only occurs to bz and 
z 3 terms. Those terms with zero phase factor are 

^±e^bz k , i p^e^z k z l z, 
l+B ' l+B 

The sum of these terms is equal to iJ k e luJkt z k = iJ k p k , where 

^) = TT§& + £M§W 2 . (3-93) 

1=1 

They are real functions bounded by Cp 2 . The other terms can be integrated. For 
example, if j G Q, j ^ k, we have 

l l^l e *^t b = df if{B) ( „ )t \ _ 

l + B 3 dt\i{u k -Uj){l + B) Fj ) 



A^ 3 )t± ( !/Cg) 

dt \i(u k + B) vyj 1 ' 



where 

error = e ^k-^)t^_ | _ v-/ ^ bp . 

and it is bounded by 

|error| < \bpj\ + \bpj\ + \bpjB\. 

Note \B\ = Cn~ 2 \b\. We can integrate other terms similarly. Summing up, we get 

~ d 

Rk = iJkPk + -j:Pk,3~ + 9kfi, (3-94) 
dt 

where p k ,3- is of the same form Q3.91| ) as R k and, adding the estimates for the inte- 
gration remainders for all five kinds of terms in ( |3.91p , we get 

\9k fi \ < n-'dbbl + b 2 \B\) + (\b\p + max \bp,\ + \bpB\) 

3 

+ ra _1 (|6|p 2 + max \bppj\ + \bp 2 B\) 

3 

+ (p 2 max \pj\ + p 3 \B\) + n~ 1 (p 3 max \pj\ + p 4 |i?|) 



3 



< p 2 max|p,| +p\b\. (3.95) 

3 
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Step 4 Final form. We now define 



qk =Pk~ Pk,i- - Pk,2- - Pk,3- 



(3.96) 



Since Pk,i- ~ nz 2 , Pk,2- ~ n 2 z 3 and Pk,3- ~ n l b 2 + bz + n 1 bz 2 + z 3 + n 1 z k z 3 , we 
have 



\Pk ~ Qk\ < \Pk,i\ + \Pk,2-\ + \Pk,3-\ < Cnp 2 . 
From (|3.89| ) and (|3.94|) we have 



q k = 22 D ki\Pi\ 2 Pk + iJkPk + Yfj=i9k,, 

1=1,..., N 



^ D kl \qi\ 2 q k + i.J k q k + 



9k, 



1=1, ... JV 



where 



(3.97) 



fl% — <7fc,l + • ■ ■ + ^,6 + <7fc,7, 

fl'fcj = Dkl (l^l 2 ^ fc _ M 2 <?fc) + iJ k{Pk - fflfc). 

i=l JV 



Since D k i = 0{n 2 ) and J k = 0(p 2 ), we have 

\9k,r\ < n 2 p 2 (np 2 ) + p 2 {np 2 ) < np 4 . 



(3.98) 
(3.99) 



(3.100) 



Collecting estimates, we have 

7 



\9k\ < yZ\9k,j\ ~ npmax\R m \ +n 2 p 2 max |pj| + [||F 5 || L i + n 2 (n 2 p 3 )] 
+ (n 2 \\r]\\ L2 + np 2 ) \F B \ + (np 4 + n 4 p 3 + n"V^) 

loc 

+ (p 2 max + p\b\) + np 4 . 



Using max m \R m \ < ||F — Fi\\ L i + (p + IMI^ )|F^| and Lemma p72| , we have 



l^fcl <np 4 + n 4 p 3 + ^11^-^11^ +(n 2 ||?7|U +np 2 )|F e | 

loc 



+ H*6|lzJ +n~ 1 pX 

loc 

< np 4 + n 4 p 3 + np(p 3 + X) + (n 2 \\ V \\ L 2 + np 2 )(p 2 + n~ l X) 

loc 

+ (np 4 + np || r/ (3) || 2 + X)+n -1 p-X 

loc 

< np 4 + n 4 p 3 + np||r/ 3) || 2 + p 2 |M| L? + n||77||* a + ||« 3 || L i 

11 11 ^l oc loc loc 11 11 -^l , 



Part 2. Ground state. 
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We have derived the main oscillatory terms of a(t) in 

a (t) = a (2) (t) + b(t), a (2) = a u z k z u 

k,ieci 

with b(t) given by ( |3.27| ). We have 



b = {c 1 Q,Im[F-F 1 }) + {c 1 Q,lmh)F e -A2,rmd, (3.101) 
with A 2 ,rmd given in ( |3.24| ), 



N , , 



iakiA-r(PkPi) 
as 



As for the excited states p^, we want to find a perturbation (3{t) of so that 
oscillatory terms no smaller than n 2 p 4 on the right side of fl3.101| ) are removed. We 
have observed in Lemma |3~2| that |6| < np 3 + nX and \(ciQ, lmh)\ < n 3 \\h\\ L i since 



h G M is almost orthogonal to Q. Hence the right side of ( |3.101| ) is of the form 

b = nizi + bj + nj l bz 2 + nj 2 b 2 z +n^r] + z^r] {2) + z 2 r] {3) + nrj 2 + (r/ 3 )io C | 
+ n 3 z Izi + b/+ n~ 1 z( J + bl) + 0(n" 2 p 4 )} + n 3 r]F e (3.102) 
+ n 2 ^ jn/ 2 + z 3 j + bzf + 0(p A /n) + nzr] + z 2 ?? + nr/ 2 + (r] 3 )i oc + ■ ■ 



Here (?7 3 )ioc means terms with same bound as ||?7 3 || L i . We shall calculate the normal 

loc 

form for b by integrating by parts those terms with orders which were crossed out. 
Notice that there are resonant terms with crossed-out orders, explicitly, terms of the 
form ri 2 |2;jfc| 2 |^| 2 . These terms cannot be integrated by parts and will remain on the 
right hand side. The final normal form equation is of the form ( p. 46 ). This procedure 



is the same as that for excited states and we shall not repeat it in details but point 
out a few key steps. 

1. There are no terms of the form b 2 , b 3 or \zk\ 2 b in the first line of ([3.102 ). Terms 
of these forms are eliminated by the Im operator. 

2. Terms of the form Cz^ziV are oscillatory if k + / ^ 0. These terms can be 



integrated. Similarly, by the observation for ( |3.50|) , terms of the form C z mi z m2 z mA1 



z m b and C z mi z m2 z m3 z m4 z ms are also oscillatory with nonzero phase factors. These 
terms can be integrated. 

3. Most terms of the form Cz mi z m2 z m:i z m4 are oscillatory. The only terms with 
zero phase factor, u mi + • ■ ■ + u m4 = 0, are of the form C | | 2 1 ^ | 2 by the observation 
for (3.50). These terms cannot be integrated and will remain in the final equation. 



Moreover, in order to survive the Im operator, these terms from F — Fi must have 
complex coefficients, i.e., they must involve t/ 2 -* and are of the order nz 2 r)( 2 \ 
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4. We need to integrate terms of the form n 2 Zjf]. They are from (ciQ, Im2AQ[(C + 
Ov + C^D i n the first line of ( ft.l02|) . Since there is only one Zj involved, the Green's 
function is approximately [Hq — eo ± (e^ — eo)] _1 and is invertible in L 2 . Hence there 
is no resonance with the continuous spectrum. This integration is carried out in 
details in Jl9, p. 193-195]. In contrast, resonant terms are of the form nz k zirf^ and 
are mentioned in point 3. 

In conclusion, we can find a perturbation f3{t) of b(t) of the form 

(3 = b + Re ^y^Cnz kl z h2 Zk 3 + Cnz k b + n 2 (z k ip k ,r])j 

+ Re 2j Cz kl z k2 z k3 z ki + Re 2J Cz k z t b (3.103) 

^fcjH \-u>k 4 ¥=0 u) k +u} t ^0 

+ Re {Cz kl ■ ■■Zkjn + Cz kl ■ ■■z k Jn 2 } , 
so that (3{t) satisfies a normal form equation, 

$= J2 B M 2 \zi\ 2 + 9b- (3-104) 

l<fc<Z<Af 

Here C denote complex constants bounded by one, and ip k denote some explicit local 
functions. B k \ are real constants bounded by n 2 , and is an error term of the form 

k 4 4 2 (H) H 2 2 2 'i 

gb ~ nz + n z +nzr] y '+nz7] + nri + nr/ + •■■, 
and we have the bounds |/3 — b\ < np 3 + n 2 p ||ry|| , 2 and 

loc 

kl < np 5 + n 4 p 4 + np 2 I|r7 (3) || r2 + n 3 p 2 \\r)\\ L2 + ™ 2 |MIl 2 + ri ll 7 7 3 |ln ■ (3-105) 

" n \ oc loc loc ^loc 

We now compute B k i, the coefficients of |^fc| 2 |^| 2 . The main contribution comes 
from 

(dQJmACV 2 )) , (3.106) 



where \( 2 r]^ is from F 5 . Although there are terms from (c\Q, lm2\^r)( 2 >) and 
(ciQ, Im2\Q{(£ + f)4 2) + ^ (2) ]) with £ = a^R + r/ 2 ), their coefficients are small of 
order 0(n A ). Since r}W = J2k,ieu U ~ X z k z iVki , and by ( ggg ) = £/*-£/* conj , 
the expression is ( 3-1061 ) is a sum of terms of the form z kl z k2 z kz z ki . Moreover, since 
rj k i has nontrivial imaginary part only if both k, I > 0, hence 



,(2) 



2j conj U\z k zir} ki - U*_z k zir} k i + irrelevant terms. 



77 

fc,i>0 fc,Z>0 



To get |zfc| 2 |^| 2 from ImA^ 2 ?/ 2 ), the relevant terms in ( 2 are v k i(z k u k ziuf) and 
^ki(z k u^ziu^), where z/ fci = 2 - d>j[.. Therefore 

(ciQ, ImACV 2) ) = ( Cl< ^' Im ^ v ki{z k ul ziuf) z k z\ conj E/^jy) 



k,l>0 



2J (ciQ, Im \v k i(z k u k ziu l ) z k Zi{-l)U*_r) k i) + (*), 



k,l>0 
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where (*) denotes terms of the form z^z^z^z^ with + ■ • • + tu^ 7^ 0. Irrelevant 
terms with + • • • + O0k 4 — are eliminated by the Im operator. Moreover, we 
can disregard the second sum since LL = 0(n 2 ) is smaller than U + . By ( 5773 ) and 
Vki = 0(n), we have 



B k i = (U+ (ciQXvkiutuj) , -Im^i) +0(n 4 

1 



A — uok — uJi — 0i 



P r A <$> 



ki 



In view of (|3.72|) - (|3.75|) , we have 



B 



ki 



-| Re^ + 0(n 4 ). 



The proof of Lemma |3.3| is complete. 



0{n A ). 



(3.107) 
Q.E.D. 



3.4 Main estimates 



Theorem |3 . 1| can be proved using the following proposition and a continuity argument. 
Recall pit) ee [p 2 + N'^onH]- 1 / 2 , A(t) = (1 + s)~ 1/2 p 2 (s), and D = 6jV|ci|7jj77 . 



Proposition 3.4 Suppose the assumptions of Theorem \3.1\ hold. Suppose for a fixed 
T > we can find the best approximation Qe(t) of ip(T), i.e., Oe(T) = in the 
decomposition (ft.l|) ofip(T) with E = E(T). Define 



(3.108) 



M t ee sup maxjpW-HELikfcWr) 172 . 2D' 1 p- 2 (t)\a(t)\, 

0<t<T ^ 

n- 4/5 P (tr 8/5 \\v(t)\\ L5 , [aw wv^wri^wb , 

loc 

2(3D)- 1 p 2 \E(T)-E \]. 

Suppose, furthermore, < 2. Then we have < 3/2. 

PROOF: Since M T < 2, we have \E{T) - £ | < 3£>p§ and, for t G [0,T], 

(Ef=il^WI 2 ) 1/2 <2pW, |a(t)|<L>p(*) 2 , 

||7 7 (t)|| L5 <2n 4 / 5 p(t) 8/5 ) 

lh (3) (*)|U <2A(t) + 2n 1/3 p 8/3 W- 

loc 

Since 77 = r/^ + rj^ 3 \ \\f]^ || 2 < np 2 and p < £o^ ? 

loc 

||r/(t)|| L? <CV + |h (3) (*)|L <C7np 2 + 2A(t). 



(3.109) 



(3.110) 



Note that A(t) < p 2 by its definition. Therefore \\r](t)\\ r 2 < Cp 2 , \\v(t)\\ r 2 nr s < 

loc loc 

Cn 2//5 p 8//5 , and the assumptions of Lemmas |3.2| and R7JI are satisfied with p = pit). 
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We have, using Lemma ^ Q and ( ^TU^j , 



\E(T) - Eo\ < \E(T) - E(0)\ + \E(0) - E \ < |[K (T) (0)| + \a Eo \] < \Dp 2 Q . 



(3.111) 



Since 



L 2 



j o\\l 2 ^ x ' 



<C 1, we have ||^(t)|| L 2 <C 1. By Holder inequality, 



\V \\ L 5/4 rZ ll'/llL2 N'/IIL 



<: llnll 2 / 3 lhll!{ 3 < o(l)(n 4 / 5 p(t) 8 / 5 ) 7 / 3 = o(l)n 28 /V 6/15 , 
< H\% 3 \\v\\f < o{l){n^p{tf/^ = o(l)n 4 V/ 3 , (3.112) 



\# < \\v\\f \\v\\% 3 < \\v\\% 3 - 4/ v /3 . 



Recall X = p 2 || 77 || i2 + n|M| 2 -2 + ||r/ 3 || L i and X < rip ||r/|| i2 + X. We have 

loc loc loc loc 

P 2 |Mb +HMli ? <P 2 (V + A)+^(V + A) 2 < np A + p 2 k + nk 2 . (3.113) 

loc loc 

Since \\rf\\ L i < p 2 |M| L 2 by fl3.112|) , we have 

loc loc 

X < np 4 + p 2 A + nA 2 , X < n 2 p 3 + npA + nA 2 . (3.114) 
We now estimate r/(t). Recall from flOg ), ( CT ) that r/(t) = U- l e~^f\{t) and 

^(t) = C/V^O) + T e - lA( *- s) P C A F„(s) ds. 

By Lemma P, (gTOj , ( gng ), ^(^(U < Cp 2 , and \\v(t)\\ L? nL s < Cn 2 / 5 p 8/5 , 



ll^ll i5 /4 < np 2 + np \\ V \\ L? + n ||r/|| 2 2 nL5 + \\ V 3 



IL5/4 



+ n 2 p 2 ||, /Nr 



< np 2 + np(p 2 ) + n(n 2 / 5 p 8/5 ) 2 + n 28 / 15 p 56/15 + n 2 p 2 (n 4 / 5 p 8/5 ) < np 2 . 

Using ( p.5| ) for 77(D) , the decay estimates of e tc and e~ itA , and the boundedness in 
Sobolev spaces of U and U~ l from Lemma |2~5|, we have 



< IKS(o) || L5 +c f\t- s|- 9/10 ||F,( S )|| L5/4 cfe 

< n 4/ 5p 8/5 + C f\t- s\~ 9 / 10 np(s) 2 ds < n 4 / 5 p(t) 8/5 + Cn^PoV 75 ^ 

for any r > 0. Here we have used p(s) = Cn~ 1 (At + s) -1 / 2 with At = CrT 2 pQ 2 and 

f |t - s|- 9 / 10 (At + s)- 1 ds < Ar r (At + t)- 9 / 10+r , (3.115) 
Jo 

for any r > 0. Taking 2r = 1/5, we get ||r/(t)|| L5 < |n 4/5 p(t) 8/5 . 

We now consider the L 2 oc estimates of rj. Recall t]^ = r]^ + • • • + r]f\ The 
estimate of nf ] is by Q and that U^e'^U = e~ ie + 0(1) [i, U] = e~ i9 + 0(n 2 ), 



V ( i\t) <(l + o(l))A(t). 
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^loc 



By the singular decay estimate ( |2.30| ) in Lemma 275 and the definition ( |3.36| ) of r] M 
with rjki = 0(n), we have 



vi 3 \t) o <^(l + t)-^«A(t), 



-3/2 



~^loc 



< 



(t — s) 3//2 n(p 2 \9\ + pmax \p k \) ds. 



By the bounds of 9 and max,,- \pj | in Lemma |3.2|, and A < p , we have 



n(p 2 |fl| + pmax |p fc |) < rip 2 (p 2 + n _1 X) + np{np 2 + X) < n 2 p 3 + npX 

k 

< n 2 p 3 + np(n 2 p 3 + npA + nA 2 ) < n 2 p 3 . 



Hence 



^loc JO 



< I (t- s)- 3/2 n 2 p 3 (s)ds<n 2 p 3 (t). 
Finally, by Lemma |3j, ( KID . \\v{t)\\ L 2 nL* ^ C™ 2/ V 8/5 , and A < p 2 , 

loc 



< p 3 + np(np 2 + A) + n^p^f + o(l)n 4 / 3 p 8 / 3 + n 2 p W /5 ) 

< p 3 + n 4/3 p 8/3 + npA. 

Hence, bounding the integrand of 774 by either L°° or L 5 -norm, we have 
V?\t) < rmin{|t- S r 3 / 2 ,|t- S r 9 / 10 }||F,, 3 ( S )|| L5/4nLl ^ 

L loc JO 

< f min {\t - sr 3/2 , \t - s|- 9 / 10 } [p 3 + n 4 / 3 p 8/3 + npA](s) rfs 
<[p 3 + n ^ 3 p 8 / 3 + npA}(t). 

Summing the estimates, we conclude 
4 

h m (t)\\ L2 <J2hf 2 ^fA + ^ + c^v^^KA^ + ^/y/ 3 )^). 



^loc 



We now estimate the error terms 5^ and Using p < e n ! ( P-48| ), ( |3.49| ), ( |3.109| ) 
( gjTTg ) and A < p 2 , we have 



\9k\ < np* + n*p 3 + np\\ri {3) \\ L2 +X 

loc 

< np 4 + n 4 p 3 + np(A + n 1/3 p 8/3 ) + (np 4 + p 2 A + nA 2 ) 



< (e 2 / 3 + n 2 )n 2 p 3 + npA. 



(3.116) 
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\g b \ < np 5 + n 4 p 4 + np 2 \\r]^\\ 2 + n 3 p 2 \\r]\\ l2 + n 2 \\ V \\ 2 l2 + n\\r] 3 \\ , 

< np 5 + n 4 p 4 + np 2 (A + n 4/3 p s/3 ) + n 3 p 2 (np 2 + A) + n 2 (np 2 + A) 2 

+ n(np 2 + A) 4 / 3 n 4 V /3 

< (el /3 + n 2 )n 2 p 4 + np 2 A. (3.117) 

We now estimate a(t). By ( |3.44| ) of Lemma |3^ we have 

\P-a\ < + |6-o| < {Cnp 3 + Cn 2 p\\r]\\ r2 ) + Cn 2 p 2 . 

loc 

Since ||r/L 2 < p 2 , we get \/3 - a\ < Cn 2 p 2 . Since a(T) = 0, we have |/3(T)| = 

loc 

|/3(T) - a(T)| < Cn 2 p 2 {T). Using flOp and ( pl7| ) for # 6 , we have 

W)\ < l/WI + / T E l*«IM 2 N 2 + \M\ds 

r T r i 

<Cn 2 p 2 (T) + / (max|5 M |)p 4 + C {e 2 Q /3 + n 2 )n 2 p 4 + np 2 A {s)ds 
Jt k ' 1 L J 

<^p 2 (t) + o(l)p 2 (t)<^p 2 (t). 

Here we have used p(s) = (p Q 2 + N~ li -y n 2 s)^ 1 / 2 and (max^j \B ki \) / (N^^^n 2 ) < D/2. 
We also have used A(s) = n 1/2 (l + s)~ 1/2 p 2 (s), p(s) = Cn -1 (At + s)" 1 / 2 with At = 
C(npo)~ 2 , and the following estimate: for £ > and At > 1, m, r > 0, m + r > 1, 

/•oo 

(Ai + s)~ m (l + s)- r ds < / (At + s)~ m ~ r ds < C(At + s)- m ~ r+1 . (3.118) 



Hence 

\a(t)\ < \(3(t)\ + Ht) - m\ < \Dp 2 {t) + Cn 2 p 2 {t) < ~L>p 2 (*). 

We now estimate the excited states Zk(t). For their initial value, we have 

\\(CE-CE )(0)\\<Cn~\p /n)\E-E \<Cn- 2 p 3 , (3.119) 

by (U) of Lemma and |# - #o| < 3.Dp§. In particular, (£* =1 ^(O) 2 ) 1 / 2 < |p 
by Let f k = \q k \ 2 . By CT) , /(0) < |pg. By O and that J fc (t) are real, we 

have 

v 

f k = J2 2ReD uhfc +2Req k9k , (k = l,...,N). (3.120) 

z=i 

Let / = fx + ■ ■ • + /n- Since | ReD^I < 57(j~n 2 , summing ( |3.120| ) over k we get 
|/| < (W^n 2 )f + 2iVmax \q k g k \ < Cn 2 p 4 + Cp((e 2 /3 + n 2 )n 2 p 3 + npA), 

k 



39 



by (JTTTBp. Let t = n~ 3 < n~ 2 p 2 . For t G [0,t ], A(t) < p§(l + 1)" 1 / 2 and 

1/(0) -/(t)|< / n 2 p A + np 2 kds < n 2 p% + np$yfc « p 2 , 
Jo 

since to "C n~ 2 p^ 2 . We conclude /(t) ~ /(0) for t G [0, to]- F° r * > ^o, we have 

£ (^f + + ^pA < (ef + n 1/2 )n 2 p 3 . (3.121) 

Since the positive part of ReD k i is bounded by n 4 and Re 25/^ < ~ 1o n2 , we have 

* 3 
f <-y j 2 l0 n 2 f 2 + Cn 4 f 2 + 2Nm a x\q k g k \<-— l0 n 2 f 2 + 2Nm a x\q k g k \. 

fe=l 

Using p 2 (t) = (pq 2 + A^- 1 7 n 2 t)- 1 and (ETI2ID , we have 
d 3 

-(p 2 ) = -iV-^onV > -— 7 o™V) 2 + 2/Vmax |g^ fc |. 
at ziV fc 

Since we have shown previously that p 2 (t) > f(t) for t G [0, t ], in particular p 2 (t ) > 

/(to), we have p 2 (t) > /(t) for all t > t by comparison principle. Thus we have 

P 2 (t) > f{t) for all t > 0. Q.E.D. 



We now prove Theorem 3.1 using Proposition 3.4 and a continuity argument. 



Proof of Theorem |37l|: By the assumption of Theorem |3J], we have M T < 3/2 
for T = 0. Denote by J the set of all T such that we can find a best approximation 
Qe(t) °f il>(T) and we have Mt < 3/2 with respect to E — E(T). Clearly J is a 
closed interval containing by the continuity of the Schrodinger equations. 

Suppose X" G J. By continuity there is a 5 > which may depend on ip and T" 
such that, for all T G [T',T' + S], there is a best approximation Qe{t) of i>(T), and 
with respect to 22 = E(T) we have Mr < 2. By Proposition |3.4| , we have Mt < 3/2. 
Hence (0, T" + 5) C J. This shows the right end of J is open and hence J = [0, oo). 

For t < T, we have |2?(t) - 2?(T)| < ||a T (t)| < Dp(t) 2 by Lemma |2j and 
Proposition [3.4| . This uniform bound shows E(t) has a unique limit 2?^ as t — > oo 
and |£(t) - E^] < Dp(t) 2 . By continuity of M T in T we have < 3/2. 

Finally we prove the lower bound ®. If (E£U IICfcH^) 172 > ip , we have 
/(0) > Cpo- We have shown previously that f(t) ~ /(0) for t < t . Summing ( |3.120|) 
over and using the error estimates, we have, for t > to, 

/ > -107o + ^ 2 / 2 + 2iVmax \q k g k \ > -12 7o + n 2 / 2 . 

k 

By a similar comparison argument we have f(t) > (/(to)~ 1 + 127 f n 2 (t — to)) -1 . Hence 
fit) > Cp 2 {t) for all t > 0. This completes the proof of Theorem |3.1| . Q.E.D. 

We remark that, since B k \ are (almost) positive, j3{t) and hence a(t) are increasing 
as t —>■ oo. This shows that the ground state gains energy from excited states even in 
the stabilization regime. 
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4 Transition Regimes 



In this section we study the dynamics in the transition regime. We will prove Theorem 



T7T| using Theorem p.l| . In this regime the natural operator is Hq = — A + V and the 



natural decomposition is 

ip = x (p H h x N (p N + £, £ G Hc(fl'o). 

4.1 Equations 

Recall H = —A + V has N + 1 simple eigenvalues e k , k = 0, . . . N, with correspond- 
ing normalized eigenvectors (p k \ H (j) k = e k cf> k , \\4> k \\ 2 = 1- We can decompose the 
solutions ip{t) = ip(t,x) respect to H : 

N 

i/>(t) = x (t) + m, x(t) = J2 ^(t)^, (4.1) 

fc=0 



where x^it) G C and G H C (H Q ). Substituting (|4.1|) into (|1.1|) , we obtain the 
following system for these components: 

ix k = e k x k + (4> k , G), (k = 0, ■ • ■ , N), 
ia£ = H Z + P C G, G = A^, 

with initial conditions x k (0) = x k and £(0) = £> From these equations we find that 
each x k (t) is oscillatory with a main oscillation factor e~ l£kt . We say that the phase 
factor of x k is — e k . Define u k {t) by 

x k (t)=e- iekt u k (t). (4.3) 

The function u k (t) has the same magnitude as x k (t) but is not as oscillatory. In 
particular, \it k \ is smaller than \x k \. We shall study the following system for u k and 



£, which is equivalent to (4.2): 



u k {t) = (0,, G(t)), (k — 0,... , N), (4.4) 

£(t) = e~ iH °% + / e - iHo{t - s) P c r 1 G(s) ds, G = \ip 2 ip. (4.5) 
Jo 

4.2 Decompositions of G and £ 

It is useful to decompose various terms according to their orders in n, so that we can 
identify the main terms. We expect that x k = 0(n) and £ = 0(n 3 ) locally after an 
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initial layer of time. We first decompose G = ip 2 ip. Using ( |4.1| ) that ip = \ + £ with 
X = J2k=o x k4 ) k, we decompose G as 

G = \ip 2 4> = X( X + 2 (X + = G 3 + G 5 + G 7 , (4.6) 

where 

G 3 = A X 2 X, (4.7) 

Gr = \ x ^ + 2X\x\X (4.8) 

G 7 = \xe + 2Xx\Z\ 2 + \et- (4.9) 

Note that G 3 = 0(n 3 ), G 5 = 0(n 5 ) and G 7 = 0(n 7 ). 

We now identify the main term of £ using the integral equation ( [4.5| ). The main 
term of the integrand is i~ l G 3 . Using ( |4.3| ) and factoring out the main phase factors 
in G 3 we have 

N 

G 3 (t)= e^-^+^um^it)^, (4.10) 

l,m,j=0 

where 

4>lmj = X4>l4>m4>j- (4.11) 

We now integrate by parts a typical term in i G 3 : 

f e -iH (t-s) p cr l ei (- ei - em+ej )s UlUmU .^ lm . ds 

Jo 

= lim rV^'f e^-^^-^u^u^P^ds 
r ^°+ Jo 



lim r l e- iHot ' 



r^o+ i(H - ei - e m + ej - ri) 

e lHos e^-tm+e^s UlUm ^ s) Pc(f)lm .] * 

We need to take a limit since H + Cj — e\ — e m may not be invertible. Define 

K{ = lim P c . (4.12) 



e im = -KLPc<Plm 3 . (4.13) 
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K\ m are bounded operators in B[L 2 ,L_ r i) with r, r' > 1/2 and r + r' > 2, see [Ffl. 
Hence we have £ 3 lm G Li r ,. We may have £/ m ^ L 2 and Im^ m 7^ only if 

e i -e z -e OT <0, (4.14) 

in particular j < Z,m. We can now rewrite the above integral as 



xix m Xj (t) xix m Xj(0) e ^ 



Irn 







Note that the choice r — > 0+, instead of r — > 0—, ensures the local decay of 
e -iH t ^ _ _ e -iH t P c imj - and of the last integral, see Lemma |2~3] . 
Summing these terms over l,m,j, we can decompose as 

= £ (2) w + £ (3) (t) = ^ (2) + (d 3) + ■ ■ • + d 3) ) , (4.i5) 

where C (^) i s th e main part of £(£), 



N 



^ 2 \t) = ■r ! r w r l (l)r L . (4.16) 

l,m,j=0 



and ^ 3 )(t) = £f } + ■ ■ • + ^ 3) is the rest, 
£ { i\t) = e- iH °%, 
d 3) (t) = -e-^ (2) (0), 



AT 



= f e- m ^- s) Per 1 (G - G 3 - A£ 2 £) 

./0 



r e -iH (ts) Pc .-i 



IIG7L1 <n|ffi +H£llK 3 IICIIK 3 , 

loc loc -^loc 

IIGs + GtIIxj <™ 2 II£IL ? , 

loc loc 

|G|| L i + max \u k \ < n 3 , 



J loc k 



\G 5 + G 7 -\et\\ LlnL5/4 <n 2 U\\ L i 



(4.17) 



We single out A|£| 2 £ since it is the only non-local term in G — G3. 
We have the following estimates for nonlinear terms. 

Lemma 4.1 Suppose, for a fixed time t, for some n < n , 

\x k (t)\ < 2n, U(t)\\ L ? n i5 < 2n, \\m\\ L * « 1. (4.18) 

loc 



(4.19) 



1 o c 
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PROOF: For the nonlocal term A£ 2 £ we have, using Holder inequality, 

||a^iu < u\\% 3 ikiiK 3 , ||A£ a eL $ u\\% 3 • (4-20) 

Similarly, ||A£ 2 £|| x <U\\% 3 U\\% 3 . Using the definition (£]| of G 7 , we have 

loc loc 

IIGVIU <nUf L? +U\\% 3 U\\% 3 - 

IOC loc "loC 



By (0) and (Q, we have ||G 3 || < n 3 and ||G 5 || < n 2 ||£|U . By 

loc 



\\G 5 + G 7 \\ L} <n*U\\ Lf +nU\\% + U\\%* U\\%* < n 2 ||£|U , 

loc loc loc loc l° c 

and||G|| Ll <||C 3 |U +||G 6 + G r |U < n 3 + n 2 < n 3 . Since |« fc | < \\G\U 

loc loc loc loc loc 

by (|4.4j) , the estimate for ii k follows. Finally, G-j — A£ 2 £ is of the form n£ 2 and 



Pr ~ A^|| LlnL5/4 <nU\\^ nL5 < n 2 U\\ LL 

loc loc loc 

by ([4.18|) . Since | \G$ H^niis/ 4 ^ ™ 2 H£Hl 2 > the ^ as ^ estimate follows. Q.E.D. 



4.3 Normal forms for equations of bound states 

Recall we write Xk(t) = e~ tekt u k (t). In this subsection we derive the normal form for 
the equations of u k , where terms of different phase factors are removed. 

Lemma 4.2 (Normal form) Suppose for some n < n Q , 

\x k (t)\ < 2n, U(t)\\ L ? ol 5 < 2n, \\t{t)\\v « l - ( 4 - 2 l) 

loc 

There are perturbations /J, k (t) of Uk{t), k = 0, 1, . . . ,N, to be defined in ( f4.39|) , satis- 
fying 

\u k (t)-fi k (t)\<Cn 3 , (4.22) 

such that 

N N 

/ifc = y^cf|//;| 2 /ifc + ^ d k ab \jj a \ 2 \n b \ 2 fi k + g k . (4.23) 

1=0 a,b=0 

Here g k are error terms, to be defined in ( |4.40| ), satisfying 

\9k(t)\<n 7 + n 2 \\^\\ L2 +nUf L? + \\t\\%' \\t\\*£ ■ (4-24) 

" ' loc loc "loc 

All coefficients cf and d k ah are bounded by one. The coefficients cf are purely imagi- 
nary. The real parts of d k ah are given by 

Re d k ab = (2 - <£)t* - 2(2 - 5 b k )j a kb , (4.25) 

where, with <fi abk = X^aMk, 

7*, = Im ( <p abk , 1— P c <j) abk ) > 0. (4.26) 

V H - e a - e b + e k - Oi J 
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Note 7^ may be positive only if e a — e& + e k < 0, in particular k < a,b. By 
assumption A2 we have 7° fc > 70 > for k > 0. The first part of Re<if 6 is positive 
and the second part negative. If k — 0, 7^ = for all a, b, hence Re(if fe > 0. If 
k = N, 7^ = for all a, b, hence Rec/f b < 0. For intermediate k the sign of ReG?f b 
depends on the relative size of these coefficients. Let f k = |/ifc| 2 , k = 0, . . . , N. Since 
^■|/i| 2 = 2Re/i/i and cf are purely imaginary, we have 

N 

f k =Yl 2 ( Rec 4) fafbfk + 2Refl k g k . (4.27) 

a,b=0 

The example 1 in §1 follows by letting N = 2, x = fio, y = fi± and z = fi2- 
Proof: Recall (O) that 



ii k = rV e *< (0 fc , G) , (k = 0, • • • , N). (4.28) 

We now proceed to derive the normal form in two steps. 
Step 1 Integration of terms of order n 3 . 

Substituting G = G 3 + G 5 + G 7 and the expression (|4.10|) for G 3 , we get 

N 

u k it)= i'W^xiXmXjWfa, <f>i mj ) + Rk (4.29) 

l,m,j=0 
N 

l,m,j=0 

R k = C l e ie <° t (<f> k , G 5 + G 7 )(t). 
By assumption A3, 

e k - ei - e m + e, ^ 0, (4.30) 
unless the two sets {k,n} and {l,m} are the same. Hence 

TV 

«*(*) = $>f N*"* + J2 i" 1 e i(efc - ei - em+,y)t u,u m %(*)(0 fc , <f> lmj ) + R k , (4.31) 

where cf are purely imaginary constants defined by 

cf = 2i~ 1 ((p k , (f) Uk ), (k^l); «4 = i -1 (0fc, 0*fcfc). 
Terms in the second group of (|4.31 ) are oscillatory and can be integrated. Define 

Uk,i = u k - u k ll 

where 



= T7 ; rUiU m Uj(t)((j) k , 4>lmj), 

tg ) <e k -e l -e m + e J ) 



45 



j^-l e i(e k -ei-e m +ej)t ^ 

9k,l(t) = - \ -— (uiU m Uj) (t)(<j) k , 4>lmj)- 

(^30|) ^ 6fc ~~ 61 ~ 6m 6j ' 

We have 

TV 

Uh,i = u k - u kl = 2Jcf \ui\ 2 u k + R k + flfjfe,i. 

We further define 

JV 

9k,2 = {\ui\ 2 u k - Ki| 2 M fc ,i} . 

1=0 

We have 

N 

u k ,i = Yl ° k i Vi + * _1 e tek *(0fc, G 5 + G 7 ) + # M + # fe)2 . (4.32) 

Step 2 Integration of terms of order n 5 . 

We next integrate 0(n 5 ) terms. We first decompose g kj i and g kt 2 since they contain 
0(n 5 ) terms. Using their explicit form we can decompose them as 

9k,i = e iekt g k>li5 + g k>1;7 , (4.33) 

9k,2 = e iekt g k , 2 , 5 + g k ,2,7, (4.34) 
where g kt ij and g k ,2,7 are higher order terms with 

Ifi'M/H < n 2 max \Rj\, \9k,2j\ < n 7 . (4.35) 
g k ,i,5 and (7^,2,5 are explicit homogeneous polynomials of degree 5 of the form 

C x h x h x h x h x h , (4.36) 

where lj G {0, . . . , iV} and C is a purely imaginary coefficient. The claim on gk,i,5 
can be seen by substituting (^4.29|) into -^(uiu m Uj) in the definition of g k:1 . The claim 
011 9k,2,5 can be seen from the definitions of g k ^ and u kl . 

We next decompose i _1 e Jefc *(0fc, G 5 ). Recall G 5 = Ax 2 £ + 2A|x| 2 £ with X = 
Y^Lo^i- Using the decomposition £ = £ (2) + £ (3) with £ (2) = J2^ m ,j=o x i x mXj £L> 
we can decompose 

C5 = (^5,1 + (^5,2 + ^5,3, 

where 

Af N 

G 5 ,i = Ax 2 ^ (-ilm£f m ) + 2A|x| 2 ^ xix m Xj (ilm^ lm ), 

l,m,j=0 !,m,j=0 
Af N 

G 5 ,2 = Ax 2 ^ xix m Xj (Ref/ m ) + 2A|x| 2 ( Re ^, 

l,m,j=0 !,m,j=0 

G 5 ,3 = Ax 2 £ (3) + 2A| X r£ (3) . 
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Recall that Im£f m 7^ only if — q — e m < 0. The term G 5i3 will be shown to 
be smaller than G51 and G^p. Both G^i and £^2 are homogeneous polynomials of 
degree 5 of the form (|4.36|) . However, the coefficients C of are real- valued Lf oc 



functions, while those of G^i are purely imaginary. 
We can now rewrite (|4.32 ) as 

N 

Uk,i = ° k i l u J,i| 2u *,i + Rk+ (i~ 1 e u, ' t (<f> k , G/5,3 + G 7 ) + g k ,i,7 + 9k,2j) , 
1=0 

where R k consist of terms of order 0(n 5 ), 

R k = e ^t {r 1 ^, G 5tl + G 5 , 2 ) + g kth5 + g k , 2}5 } . (4.37) 

Inside the bracket are polynomials in xi and Xi, I — 0, . . . , iV. The only terms with real 
coefficients are those from G^\. As in Step 1, the above expression can be separated 
into two groups: those with zero phase factor and those with non-zero phase factors. 
The phase factor of a typical term in the above expression is of the form 

e,t - e h - e h - e h + e h + e h , 

where lj G {0, . . . , iV}. To get zero, we must have the form e k — e k — e a — e b + e a + e b 
for some a, b by assumption A3. Thus a typical term in the first group is of the form 
Ce iekt x k x a x b x a x b = C\u a \ 2 \u b \ 2 u k . 

We can integrate by parts those terms with non-zero phase factors as in Step 1. 
The other terms remain. Hence we can rewrite R k as 

~ N d 

Rk = d ab \Ua\ 2 \u b \ 2 U k + -Tt(%, 2 ) + 9k,3, 



a,b=0 

]k 



for some order one constants d^ b , some explicit homogeneous polynomials u k 2 of 
degree 5, and remainder terms g k $ with 

\g kfi \ < n 4 max \uj\. (4.38) 



3 

]k 



We will identify RecQ, in a moment. Define 

ilk = u k ,i - u k 2 = u k - uZ x - u k 2 . (4.39) 



We have 



[ik = u kA - u k 2 

N N 
= 2jcf |Mi,l| 2 Mfe,l + 2J ^ab \u a \ 2 \u b \ 2 U k + 9k,3 
1=0 a,b=0 

+ (r l e iekt ((j) k , G 5j3 + G 7 ) + g kX7 + g kt2J ) 

N N 

= ^2c k \fll\ 2 fl k + Y d ab \^a\ 2 \Vb\ 2 Vk+ 9k, 
1=0 a,b=0 
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where 



g k = i 1 e iekt {(f) k , G 5 ,3 + G 7 ) + g k ,i,7 + g k ,2j + 9k,3 + 9k^ (4.40) 



N N 



9k,4 = {\ll ltl \ 2 U ktl - \Hl\ 2 fl k } + ^ d ab {\Ua\ 2 \Ub\ 2 Uk ~ \f^a\ 2 \t*b\ V* } • 

1=0 a,b=0 

We have \gk,s\ i$ n "' '■ Collecting estimates, using the explicit form of £75,3, \R k \ < 



11^5 + G7 1| l 1 ' anc ^ Lemma JO, we get 



loc 



l#fc| ~ 11^5,311 ri + IIG7L1 +n max\Rj\+n +n max\uj\+n 

loc loc j j 

<n 7 + n 2 ||e (3) || L2 +n*U\\ Lf +n\\£\\l> + I^L ■ 

loc l° c loc loc 



Using ||^|| L 2 < n 3 + ||^ 3 ^|L 2 > we can remove n 4 ||£||r2 and get ( |4.24| ). 

lor 11 ' ' ^^^^ Inn 



We have noted that cf are purely imaginary. We now compute the real parts of 
d^ b . From the previous discussions on gk,i,5, g k ,2,5, G^,i and £5,2, we know that the 
only contribution to Rec^ b is from i~ l e iekt {(j) kl G51). We have 



/ ''<"'■'(< >/.,, G 5> i) = 2J -e l£fc * x a x b xix m Xj A0 a 6 Im f/ m ) 



.Vfc) ^5,1; — 



+ 2j 2ei6fc< XaXbXlXmXj (0fc, A0 a 6 1H1 £/ m ) , 

and we are interested in those terms of the form Ce iekt \u a \ 2 \u b \ 2 u k . As mentioned, 
Im£j? m 7^ only if ej — e\ — e m < 0, in particular j ^ l,m (not sufficient though). 
Therefore, to get products in the first summation we need j = k and the two sets 
{/, m} and {a, b} are the same; in the second summation we need j = a and the two 
sets {l,m} and {k,b} are the same. Hence ^2ab=o0^ e dab) \ u a\ 2 \ub\ 2 Uk is equal to 

N 

"(2 - S b a ) \x a \ 2 \x b \ 2 u k (0,, A0 a 6 lm£ a \) 

a,b=0 
N 



2(2 - 5 b ) \x a \ 2 \x b \ 2 u k (0 fc , A0 a fe Im4 a 6 ). 



a,b=0 



Note we have the factors (2 — 8 a ) and (2 — <5y since we have two choices for assigning 
I and m if a ^ b (or k ^ b). Recall \<p a <Pb<Pk = (fiabk, CL = - K im^4>imj, and 
tL = (0/mjJmif/ m P c <^ mi ). We conclude 



Re<4 = -(2 - 8 b a )(<f> abk , Im&) + 2(2 - 5 b k )(<f> kba , Im&) 

Q.E.D. 



(2 - - 2(2 - S b ) 7 a k 



48 



4.4 Main estimates 



In this subsection we prove estimates in the transition regime for the solution ip(t). 
The following proposition is the main result of this section. It proves the same 
result of Theorem |1.1| under weaker assumptions (|4.41|) - (|4.42 ) below. Note that the 
assumption ||£o||y < w/2 in ( |1.9|) of Theorem |1.1| implies, using Lemma |2.3| , 



-itH 



-9/10 



-itH 



Zo\\ L 2 <n(l + t) 

loc 



-3/2 



Hence it implies ( 4.42|) . Eq. ( 4.42|) is motivated by ( |4.46|) and P2] , and may be more 
convenient for future application when we study the asymptotic profiles of all small 
solutions. See for two bound states case. Also note that, by ( 4.41 ) below, we have 
2 = HV'ollz^ <C 1 and hence max^ 11^(^)11^2 -C 1. However, we do not 



assume any bound of ||^o|Il2 in terms of n. This is essential for future application. 

Notice that we define t 2 by ( [4.45| ) to be the time when the size of the excited 
states decays to the order n 1+a . This is more than sufficient since we only need it to 



be smaller than Equ in order to apply Theorem I3J. 



Proposition 4.3 Assume the assumptions A0-A3 given in §1. There is a small 
constant n > such that the following hold. Suppose ip(t) = ip(t,x) is a solution of 
( |1.1|) with the initial data ip(0) = ipo = XqVo + • • • + x%(pN + £o satisfying, for some 
n < n , < a < 1/10, 



| n < \x, 



3 n 2 

Xq\ > 77 



1 2 



„0 12 



01 

3-e 



< n 2 



1 < 1, 



\x 



1 2 



1/2 



> 2\x 



Oil+o- 

ol > 



+ • • • + \%]\[\ 

and for some t* G [1, n~ 4 ~ 2cr ], for all t > 0, 

lle-^olL < A 5 (t) = n(l + t)- 9/10 + Cn%(U + t)~ 9/1 °, 



(4.41) 



\e~ UHo Zo\\ L 2 < A(t) = n(l + t)- 3/2 + Cn 3 



Let to = n 2 . There exist t 2 >t\> t such that 



U + t 



[l + t) 



-1/2 



(4.42) 



ti < C 4 n- 4 log -, C^n" 4 <t 2 < C 4 n 
n 



-4-2ct 



(4.43) 



for some constant C4 > 1 independent of n, and we have, for < t < t 2 , 

\x (t)\ > f sup \x (s)\, 



0<s<t 



\x k (t)\<in, (k = 0,...,N), 
U(t)\\ L s<C 5 nh^ + ^A 5 (t), 



(4.44) 



urn 
\e 3) (t)\\ 



L 2 < C 5 n 3 



1 2 



<C 5 n 5 } 



A(t), 

for t > t , 
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where C5 is an explicit constant. We also have, att — t\ and t = t 2 , 

l^o(*i)| a , ko(* a )| a > \l- 2 - N n\ (Ef=il^fe)| 2 ) 1/2 < Mt 2 )| 1+<J . (4-45) 
Moreover, we have the following out-going estimates for £(£2)-' For T > 0, 

m\\ L5 <Cn%U+t 2 )(t 2 + T)-^ 10 , 

^°£(t 2 )|L 2 <Cn 3 ^ (1 + rH/ 2 (U(>) 

^loc to + T 



-i-ri/b < 



PROOF: We will prove estimates ( |4.43|) -( [4.45| ) using a continuity argument. Hence 
we can assume the following weaker estimates: For < t < t 2 : 



|x (t)| > \ sup |x (s)|, 

0<s<t 



\x k {t)\<\n, (k = 0,...,N), 
||£(t)|| L5 <2C 5 nV/ 10 + 2A 5 (t), (4.47) 
||£(t)|U < 2C 5 n 3 + 2A(t), 

loc 

||£ (3) W|| L2 <2C 5 n 5 , fort>t . 



If we can prove (|4.43| )- (|4T45|) assuming ( [4.471 ), then ([4.43|) -( ^45|) hold true for all 



t G [0,^] since (|4.47|) is never violated. By definition we have As(t),A(t) < n. By 
( |4.47| ) we have ||£(£)|| L 2 nL5 < 2n. Hence the assumptions of Lemmas [4. 1| and Lemma 

loc 

fO| are satisfied. Since ||?/'(t)|| i 2 = HV'oll^ <C 1, we have ||£(£)|| L 2 <S 1- Also, since 
t<t 2 < C 4 n'^- 2a , 

UWifi Z n % /W + "(1 + ty 9/W + CnH\l w < n 13 ^~ a / 5 + n{\ + t)" 9 / 10 , (4.48) 
||£|| L 2 < n 3 + n(l + t)~ 3/2 + Cn 3 (l + t)' 1 / 2 < n 3 + n(l + t)^ 2 . (4.49) 

loc 

We now estimate £(£). Recall (|4.5| ). 

f (t) = e~ itHo Co + / e- i( *- s)Ho P c r 1 G(s) ds. 
Jo 

Note ||G|| L5/4 < ||G3|| L s/4 + ||G - G? s - A£ 2 £|| L5/4 + ||A£ 2 £|| L5/4 and ||G 3 || L5/4 < Cn 3 . 
Using Lemma [O], (|4.48| ) and (|4.49|) , 



G -G 3 - A£ 2 £|L lnL5/4 < n 2 ||£|| L 2 < n 2 (n 3 + n(l + S )" 3 / 2 ), (4.50) 



2 ^"x,5/4 < neii^iieiiK 3 <o(i)[^ 13/5 - ff/5 + n(i + s )- 9 / 10 ] 7/3 - (4.5i) 
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Hence \\G{s) || i5/4 < n 3 + n 7 ' 3 (s)" 21/10 and 



lie(*)ll L B < ||e-^°6|| L5 +^ / \t-sr'"\\G(s)\\ L5/i ds 



< A 5 (t) + C / |t - s|- 9 / 10 [n 3 + n 7 / 3 < f A 5 (t) + CnH 1 / 10 . 

Jo 

Hence ||£(£)|| L 5 is estimated. For the L 2 oc norm, we use the decomposition ( [4.1 5|) that 

£( t )=£(2) + £(3) with ^(3) 

By ( WM) for £ we have 

if\t) ^ <A(f). 



£(t) = ^ +£® with ^ = J2 5 n=! d 3) - We have ||£ (2) |L < Cn 3 by its explicit form. 

J J Mi/, 



^loc 



By the singular decay estimate (|2.5|) in Lemma |2.3| and the definition ( |4.13| ) of £; 

we have 

$\t) 2 < Cn 3 (l + t)- 3 / 2 < Cn 2 A(t). 

"^loc 



Using ( |2.5| ), ( 4.13 ) again and the bound of max.,- \v,j\ in Lemma 4~T, we have 



. <C/ (*-s)~^n 2 max|u 3 -|ds < C / (t-s)-" /z n 5 ds < Cn 5 . 



^loc 



-3/2 2 



"3/2^5, 



For ^^(i), bounding its integrand by either L°° or L 5 -norm and using (|4.50| ), we have 



< 



< 



f min { |t - s\-*'\ \t - S |- 9 / 10 } ||(G - G 3 - \eO(s)\\ LlnL5/4 
Jo 

f min { \t - S |- 3/2 , \t - S |- 9/10 } n 2 (n 3 + n{\ + s)- 3/2 ) 
Jo 

< n 5 + n 2 A(t). 

Finally we give two estimates for £5 . By Holder inequality and (|4.48|) , 

\\e\u L5/ , < u\\V < 0(1)^-^ + n (i + ,)- 9 / 10 ) 5 / 3 . 

< n 13/3- ff /3 +ri 5/3 (l+s) -3/2_ 

The same estimate for ^ shows 

Zs\t) 2 <n 13/3 ^ /3 +n 2/3 A(£). 

"^loc 



ds 



(4.52) 



Summing the estimates we conclude 



loc IOC ' II f- 2 



J = l 



; 2 
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(3) 

We now estimate £5 in another way. We have 



\m\ LW/ . < llellX 9 nL5 U\\T < o(l)(n^-^ + n(l + s) -^ )2 o/9 

< n 52/9-4 CT /9 +n 20/9 (1 + s) ^ 

Bounding the integrand of ^ (t) by either L 6 or L 5 norm we have 



^loc 



< 



< 



;s)l lL 6 / 5 nL 5 /4 ds 



[ min{|t- S |-\ |t- s |- 9 / 10 }||A^ 
Jo 

{\t-s\-\ \t-s\- 9 / w }(n^ 9 -^ 9 +n 2 °/ 9 (l + sr 2 ) ds 



mm 



< n 52/9~4 CT /9 lQg(1 + h) + n 20/9 (l + ^-1 j og(1 + ^ 



Summing the estimates and using ti < Ca£ 



-4-2o- 



we get 



ie (3) IL L <E|kf f2 <Cn 5 + |A(t)+n 2 °/ 9 (l+t)- 1 log(l + t 2 ). 



"^loc 



Clearly we have ||^ 3) (t) II f2 < Cn 5 for t > t 



o = n" 



We now estimate the error term gk{t) in the equation ( [4.23|) of Uk- Recall (|4.24j ), 

|^(t)|<n 7 + n 2 ||^|| L2 + n\\t\\% +U\\% 3 UWf- 



Using (gg7|)-(ggg), A 5 (t), A(t) < n, and the computation in (fO)g ), 

l<7fc(*)l < + ||^ (3) || L2 + ^(n 3 + A) 2 + (n 3 + A)V3 ^13/3-^/3 + n 2/3 A(s)) 

loc 

<n 7 + n 2 \\^ ] \\ L2 + n 14/3 " CT/3 A + nA 2 . 



For t < to = n 2 we use the first estimate of 



(3) 



V 5 



(3) 



and we have ||£^|| r2 < 

L? MM L 1qc 

loc 



< n i3/3-<j/3 _|_ ^ Hence 

\9k(t)\ <n 6+1 ^ 3 + n 2 A, (te[0,t o })- 
For t > t we have ||£ (3) 0)|| L 2 < Cn 5 and A(t) < n 2 . Thus 

loc 

I^WI<n 7 , (te[t ,t 2 }). 



(4.53) 



(4.54) 



We now estimate bound states. As illustrated in Example 1 of §1, some bound 
states may grow at intermediate time. We will use the following monotonicity prop- 
erties. Recall that f k = \fik\ 2 , k = 0,... ,N, satisfy ( |4.27| ). Motivated by ( |1.19| ), we 
consider the following functions: 

N N 

/(t) = /o(t), g(t) = J2fk, g+(t) = J2 2 ~ k h- ( 455 ) 



k=l 



k=l 
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By ( TO) , ( |Og ) and (^1) we have 

/(0) > 2-V- 2e , 2n 2 >/(0) + ( 7 (0)>2- 1 n 2 , /(0) + <7+(0) > 2- 1 ~ A V, (4.56) 

if n is sufficiently small. Although g and g + are comparable, (/ + g)(t) is almost 
monotone decreasing while (/ + g + )(t) is almost monotone increasing in the following 
sense: 

^(/ + < 2(N + f ) max |^ fc |, 
at k 

-(f + g + )(t) > -2{N + l)max\ji k g k \. 
dt k 

We now prove ( ^57|) . By flOED and fl£27D , 

N N 



(4.57) 



^(/ + ^)-^2Re/i fc ^= E 2(Re d k ab )f a f b f k 

k=0 a,b,k=0 

N 

= 2 K 2 - <M 6 " 2(2 - <S*)t&] /Jtf*. (4.58) 

a,b,k=0 

Switching a and k in the terms with minus sign, we see that the above sum is non- 
positive. Similarly, 

, N N 

Jt (f + g+)-J2 2l ' kRe ^9k= E ^- k (Red k ab )f a f b f k 

k=0 a,b,k=0 

N 



E ^ K 2 - ^ - 2 (2 " ^* 



a,fe,fc=0 
iV 



= E [ 2 ^( 2 -M-2 2 - a (2-M]/ a / 6 / fc . 

a,b,k=0 

In the last line we have switched a and k for terms with minus sign. If k > a then 
= 0. Hence 2 1 ~ k > 2 2 ~ a for nonzero terms and the above sum is non-negative. 

We have shown (|4.57|) . We now estimate the bound states in three steps. 

Step 1. Initial layer. 

In this period the dispersive part disperses away so much that it becomes negligible 

locally. The time it takes is less than i — n~ 2 , which is not long enough to change 

the magnitudes of the bound states. Explicitly, by (|4.27|) , ( f4.47|) and ( |4.53| ) we have 

\f k (s)\ < Cn 4 a 2 + a ( n 6+1 / 3 - ff / 3 + n 2 A(s)) , 



for k = 0, 1, . . . ,N, where a = \J f k (0). Thus, for t e [0, t ] 



\fk(t) - f k (0)\ 



< 



[ nW + a (n 6+1 / 3 -^ 3 + n 2 A(s)) ds 
Jo 

< n A a 2 t + an Q+1/3 - a/ H + an 2 (n + n 3 Vu) 

< n 2 a 2 + n 3 a < n e a 2 + n 6_e . 
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Since f k (0) > n 6 ~ 2 % we get \f k (t) - f k (0)\ « / fc (0). 

Step 2. Transition regime (i). 

After time t = t most mass of the dispersive wave is far away and has no effect 
on the local dynamics. Two possible situations can occur: We either have /(to) > 
2~ 2 ~ N n 2 or the opposite. In the first case we define t\ = t and jump to next step. 
We now focus on the second case. In the second case the ground state begins to grow 
exponentially until it is of order n. The time it takes is of order n~ 4 logn _1 . Define 

t x = inf {t : f (t) > 2~ 2 - N n 2 } . (4.59) 

t>to 

By assumption t± > t . We want to show that 

to < h < t[ =t + 2 8+2N N- f0 - 1 n' 4 log - . (4.60) 

n 



Suppose QOOP fails, that is, f (t) < 2- 2 ~ N n 2 for all t < t[. We have by ( [£55D , ( ^54|) 



and (|C57D , for all t G [t ,*'i] 



9+(t) > (f + g+)(t ) - f Q (t) - / max | H^lds 

J to 

> 2 - 1 - N n 2 - 2- 2 - N n 2 - Cnn 7 > 2^- N n 2 . 



In particular, the coefficient of the linear term in ( [4.27 ) of /o has a lower bound 



N N N 



2(Red° 6 )/ a / 6 = 2 ( 2 " 5 a)7 a °Ja/ fe > ^2 7o / a 2 > 2 7o ArV (4.61) 

a,b=0 a, 6=0 a=l 

> 2 7o ArV > 2- 5 - 2 ^ 7o iV- 1 n 4 . 



Since y/fitf) = /o(*)/M*)| < 3/ (*)/|a*o(0)| < 4/ (tK 3+e by the first line of O, 
we have for t >t , 

f (t) > 2~ 5 - 2 %iV- 1 nVo - V^^ 7 > 2~ 6 - 2JV 7o iV- 1 n 4 /o, 
if n is sufficiently small. Hence 

/o(t) > /o(*o) exp {2- 6 ~ 2 %iV- 1 n 4 (t - t )} , (4.62) 
for t £ [to,ti]- However, using the definition of t[ and 

/o(to)>!ko(to)| 2 >^ko(0)| 2 >^n 6 - 2e , 

( |1.62j ) implies Jo(^i) > 2~ 2 ~ N n 2 , which is a contradiction. This contradiction shows 
the existence of t x E [t ,t[] so that / (ti) = 2~ 2_A n 2 . We also have and ( P%D 

for all t <t\. 
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Step 3. Transition regime (ii). 

For t > ti, fo(t) is large enough to cause g(t) to decay in a rate we can control. 
Define 

t 2 = inf {t : g(t) < n 2+2a \ . (4.63) 

We want to show that 

h < t 2 < t' 2 = ti + n" 4 " 2 * 7 . (4.64) 

Suppose the contrary, then git) > n 2+2cr for all t < t' 2 . By (^4.47|) , for t > t\ 

we have ^Mf) < /o(^)(|/o(ti)) _1/2 < 2^ 3+N ^ 2 nf (t). Recall from flOID that 
£f 6 =o 2 (Re > 2 7o iV-^ 2 . Hence for t G 



2j> 



/o(*) > 2 l0 N- L g 2 - V/oCn 7 > 0. 



Thus f > and / (*) > 2~ 2 - A n 2 for t G In particular (fO§ holds for t > £ a 

and ( |4.45| ) holds for Xofe) should t 2 exist. By eq. (|4.27| ) for f , 



g(t) < ~fo(t) + Ef=o 2Re Vk9k 

= ~ {Ea, b =o 2 (^d° ab )f a f b } / + Ef=i2Re/2^ fe 

< -2 7o iV- 1 ( 7 2 /o + Ef=i2Re^ fe 

< -2- l - N n 2 l0 N- l g 2 + CW 7 

< -2- 2 -%iV" W, 

if n is sufficiently small. In the second line we have cancelled 2 Re/io^o- I n the third 
line we used ( [4.611 ). Thus, by comparison principle, 

g(t) < [9(h)- 1 + 2- 2 -%iV- V(t - h)}- 1 , (h<t< 4). 

for t G (ti,t' 2 ]- However, by the definition of t 2 this implies g(t' 2 ) < n 2+2cr , which 
is a contradiction. This contradiction shows the existence of t 2 satisfying ( |4.64j ). 
Similarly, since g(t) > —Cn 2 g(t) 2 for t\ < t < t 2 , we have a lower bound g(t) > 
[<7(*i) _1 + Cn 2 [t — t^]^ 1 by comparison principle. This implies t 2 — t\ > Cn~ 4 . We 
have proven ( |4.43| )- ([4.45|) in Proposition [O] using assumption ([4.47|) . Since ( |4.47| ) 



holds for t = 0, it holds for all t < t 2 by continuity. Hence ( |4.43| )- (|4.45| ) are proven. 



Finally we prove ( [4.46[ ) for e tTH °£,(t 2 ). Let t = t 2 + r. We have 



°£o + J{t), J(t) = f 2 e -*('-'>*° P c G(s) ds. 

Jo 
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The estimate of e~ itH °^ is by (|Og ), since A 5 (t) < CVi 3 (i* + t 2 )(l + i)~ 9/10 and 
A(i) < <7n 3 (t* + t 2 )(l + + r)- 1 / 2 if t 2 > n~ 2 . To estimate J(t), we use the 

following integral inequalities: For t > T > 1, 

[ T \t-s\- 9 / 10 ds<CTt- 9/1 °. (4.65) 

[ T mm{(t - S )~ 3/2 , (t - s)- 9/10 } rfs < CTr 1 (t - T>~ 1/2 . (4.66) 
Jo 



See p2|, Lemma 2.6]. Since ||G(s)|| L i nL5 /4 < Cn 3 for s G [0,t 2 ], we have 

\\J(t)\\ L5 < C f 2 \t - s\~ 9 / 10 \\G^s)\\ L5/4 ds<C f 2 \t - s\- 9 / 10 n 3 ds 
Jo Jo 

<Cn%{l + t)- 9 / 10 . 

\\J(t)\\ Lf <C T min{\t-s\- 3 '\ |t- S |- 9/10 }||G( S )|| LlnL5/4 da 
Jo 

< C [ 2 mm{\t- s\- 3/2 , \t- s\- 9/10 }n 3 ds 

Jo 

< Cnhzt^tl+t-tz)- 112 . 

The proof of Proposition [O] is complete. Q.E.D. 



Proof of Theorem |1.1| To prove Theorem |T7T] using Theorem [3.1| , we want 



to use ijj{t 2 ) as initial data and decompose ip(t 2 ) according to (|3.1|) , i.e., with respect 
to C, and check the assumptions (|3.4|) - (|3.5|) . If the last condition in (|4.41 ) fails, we 
simply set t 2 = 0. Let p = {\xi{t 2 )\ 2 + h |xAr(t 2 )| 2 ) 1 ^ 2 - We have 

M* 2 )|~n, p = (l + o(l)K +ff , ||£(* 2 )|U <Cn 3 . (4.67) 

loc 

We may redefine n = \x (t 2 ) |. Let Qe be the unique nonlinear ground state so that 
Qe = W>o + k with k _L 0o, k = 0{n 3 ). Also choose e te ° = n^xofo). We have 
\\ip(h) ~ QE e i9o \\ r2 < \\k\\ L 2 + \\ip(h) - x (t 2 )4> \\ L 2 < n 3 + p + n 3 < p , which 
is much smaller than n. By Lemma [2.2| , for each close to Eq, we can decompose 
ip{t 2 ) as VK^) = [<5_b + uRe + C + v) e%d using the decomposition ( J3 . 1| ) with respect 
to We have 

e^fy^) -Q E = aR + C + 7] = k + Z + ^ (4.68) 

where Z = e-*[ari(t 2 )0i + • • • + a;jv(t 2 )] and £ = e^^ta)- With £ = £ in ( ggg ), we 
have a^ = (ciQ, Re A; + Z + £) = Ci(n0 o + A;, Re /c + Z + £). Since O -L + Z + £), 
we have 

\a Eo \ < n 3 (n 3 + p + ™ 3 ) < Cn 3 p , 
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which is much smaller than Pq. Denote Pe = J2k=i^k- We also have 
||C - PeZ\\ = \\P E (k + OH < n 3 + n 2 U\\ Lf < n\ 

loc 

and hence ||£|| = Po + 0(n 3 ). We have verified (|3.4j) with P = P . F° r an P with 
|P - P | < £>Po, write Vi = Pc C (k + Z) and rj 2 = P c c £ = P c c e~ ie £(t 2 ) . We have 

e sC r] = e sC rn + e sC rj 2 . 

Since rji is a local function bounded by \\k\\ +n 2 \\Z\\ < n 3 , we get 

||e s SiL s < Cn 3 (l + S )- 9 / 10 « n^p{sf\ 
\\e sC Vi\\ r i < Cn 3 (l + s)- 3 / 2 « n l ' 2 {l + sy l / 2 p 2 {s). 

11 11 ^loc 

Recall p(s) = [pq 2 + N~ l -f n 2 s}- 1 ^ 2 with p ~ ™ 1+fT - Write £ as £ = -i(H Q -E)+W, 
where Wh = 2XQ 2 h + XQ 2 h is a local operator bounded by n 2 . We have 

eS? 2 = P c c e^ H °- E h- i6 i(t 2 ) + f ' e {t ~ s)c We- is{H °- E) e- w £(t 2 )ds. 

Jo 

Hence 

\\e t %\\ LB <\\e- itH °t(t 2 )\\ L5 + f\t-s\-^n 2 \\e- H ^{t2)\\ L , ds, 

Jo loc 



e tC V2\\ L2 <\\e- itHo £(t 

1 1 -^i^^ 1 1 



2j 



+ fmmilt-sr 3 / 2 , |t- s |- 9 / 10 }n 2 ||e-^(t 2 )| ly 



Using the out-going estimates Q4.46|) for e tsH °£(t 2 ), we conclude the estimates (|3.5|) 



for for all P close to E with |P - P | < -DPo- Q.E.D. 
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